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ABSTRACT 


The combinatorial geometry allows efficient Monte Carlo 
radiation transport calculations of detailed three-dimen¬ 
sional geometries. Because regions and media are formed 
by combination of basic bodies such as boxes, spheres, 
cylinders and others, the input required of the user is both 
relatively simple and easily modified. The M0RSE code is 
a multigroup neutron and gamma ray transport Monte 
Carlo code that may solve either neutron, gamma ray, or 
coupled neutron-gamma ray problems in either the forward 
or adjoint mode. M0RSE has a wide variety of available 
input options, including splitting, Russian roulette, expo¬ 
nential transform, energy biasing, importance regions, 
albedo surfaces, and the scoring options available in the 
SAMB0 analysis package. 

This document details the incorporation of a new version 
of combinatorial geometry into the M0RSE code and is 
meant as a user's manual. All new subroutines, including 
the combinatorial geometry, are documented in detail, in¬ 
cluding flow charts. Input instructions for all M0RSE 
options and layouts of all M0RSE commons are given for 
user reference. The PICTURE code has been modified to 
work with the combinatorial geometry, allowing two-dimen¬ 
sional printer pictures of combinatorial geometry models. 
Two M0RSE sample problems and a PICTURE sample 
problem are discussed. Input and output from these prob¬ 
lems and the complete code are available from the 
Radiation Shielding Information Center. 
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I. INTRODUCTION « 

(i) 

The M0RSE Monte Carlo code was first distributed'*' in 

I (2 2 ) 

early 1970 with a revised versions ' ’ released in the fall of 1970. 

There have been several additional features^’ added in the past year. 

\ 

This report describes the addition of a version of the MAGI developed 
combinatorial geometry' ’ (CG) to the M0RSE code and also serves to 
document a version which is operational on the UNIVAC 1108. The 
definition of variables in labelled commons, layout of blank common 
and input data have been collected for easy reference. A description 
of the PICTURE program for the combinatorial geometry package and 
descriptions of routines required to solve two sample problems are 
included. , 

l 

The replacement of the 05R geometry module with the CG 
module required significant modification of the CG package. These 
modifications were required so that the CG ipodule could be added to 
M0RSE rather than adding MORSE to the CG package. Only those 
M0RSE routines that have an interface with the geometry package were 
changed and all options in M0RSE are available with the CG module. 

To retain all M0RSE options it wa^ necessary to write several new 
routines as part of the CG module. 

V Features of M0RSE which are dependent on the geometry 

module include the ability to (1) deter mine, both the media and an im¬ 
portance region given a set of coordinates, (2) to track a particle through 
the system including internal voids and surrounded by an external void, 
i (3) determine the number of mean free paths between any two sets of co¬ 
ordinates (one may be in either an internal or an external void), (4) 
determine the vector normal to a surface at any point, and (5) permit 
particle reflection at a boundary without "getting lost. ” Table I lists 
the M0RSE module routines which were modified to permit an interface 
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TABLE I. 


M0HSE Routines, that Depend on Which 
Geometry Module is Used 


EUCLID (3) 
G0MST (2) 
INPUT 1 
MS0UR (1) 


with the! CG module. The numbers after the routines Relate the features 
mentioned above with the function of the routines. The routines without 
numbers are concerned with the input of geometry information. Table n 

TABLE II. 


Geometry Module Routines 


05R Spherical Geometry 

Combinatorial Geometry 

GE0M 

(2,3) 

GENI 


G0MFLP 

(5) 

GTVLIN 


J0MIN 


GG 

(2,3) 

'j L00KZ 

(1) 

G1 

(2,3) 

N0RML 

GTV0L 

(4) 

ALBERT 

PR 




G0MFLP 

(5) 


' 

J0MIN 



\ 

i 

L00KZ 

(D 


N0RML 

GTV0L 

(4) 
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lists the geometry module routines for the 05R spherical and the com¬ 
binatorial geometry packages. 

Only minor changes were made to INPUT and MS0UR and no 
additional documentation is provided. INPUT was split into two routines, 
INPUT1 and INPUT2, to allow a more efficient overlay. Numerous 
changes were required for subroutine EUCLID and G0MST and new de¬ 
scriptions and flow charts of these routines aye given. A complete 
description of the combinatorial geometry package is included and 
sample problems are provided. 
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II. COMBINATORIAL GEOMETRY 

2.1 Body Types 

Combinatorial geometry (CG) describes general three dimen¬ 
sional material configurations by considering unions, differences 
intersections of simple bodies such as spheres, boxes, cylinders, etc. 

In effect, the geometric description subdivides the problem space into 
* 

unique zones. Each zone is the result of combining one or more of 
the following geometric bodies. 

1. Rectangular Parallelepiped (RPP) 

2. Box (An RPP randomly oriented in space) 

3. Sphere 

4. Right Circular Cylinder 

5. Right Elliptic Cylinder 

6. Truncated Right Angle Cone 

7. Ellipsoid 

8. Right Angle Wedge 

9. Arbitrary Convex Polyhedron of 4, 5, or 6 sides. 

Body types S-9 may be arbitrarily oriented with respect to 
the x, y, z coordinate axes used to determine the space. Body 1, a 
special body described below, must have sides which are parallel to 
the coordinate axes. 

The basic technique for the description of the geometry con¬ 
sists of defining the location and shape of the various zones in terms of 


To avoid confusion between importance regions and combinatorial 
geometry regions, we depart from previous combinatorial geo¬ 
metry descriptions and use the term zone to indicate a combinatorial 
geometry region which is designated by the variable IR. The term 
region is reserved for an importance region. Thus the zone index 
is IR. 
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the intersections and unions of the geometric bodies. A special oper¬ 
ator notation involving the symbols (+), (-), and (0R) is used to describe 
the intersections and unions. These symbols are used by the program 
to construct information relating material descriptions to the body de¬ 
finitions. 

If a body appears in a zone description with a (+) operator, it 
means that the zone being described is wholly contained in the body. If 
a body appears in a zone description with a (-) operator, it means that 
the zone being described is wholly outside the body. If the body appears 
with an (0R) operator, it means that the zone being described includes 
all points in the body. In some instances, a zone may be described in 
terms of subzones lumped together by (0R) statements. When (0R) 
operators are used there are always two or more of them, and they re¬ 
fer to all body numbers following them, either (f) or (-). 

Techniques for describing a particular geometry are best il¬ 
lustrated by examples. Consider an object composed of a sphere and a 
cylinder as shown in Fig. 1. To describe the object, one takes a spher¬ 
ical body (2) penetrated by a cylindrical body (3) (see Fig. 1). If the 
materials in the sphere and cylinder are the same, then they can be 
considered as one zone, say zone I (Fig. lc). The description of zone I 
would be 

I = 0R + 2 0R + 3 

This means that a point is in zone I if it is either inside body 2 or inside 
body 3. 











If different materials are used in the sphere and cylinder, 
then the sphere with a cylindrical hole in it would be given a different 
zone number (say J) from that of the cylinder (K). 

The description of zone J would be (Fig. Id): 

J = + 2 - 3 

This means that points in zone J are all those points inside body 2 
which are not inside body 3. 

The description of zone K is simply (Fig. le): 

K = + 3 


That is, all points in zone K lie inside body 3. 

Combinations of more than two bodies and similar zone de¬ 
scriptions could contain a long string of (+), (-), and (0R) operators. 

It is important however to remember that every spatial point in the 
geometry must be located in one and only one zone. 

As a more complicated example of the use of the (0R) operator, 
consider the system shown in Fig. 2 consisting of the shaded zone A and 
the unshaded zone B. These zones can be described by the two B0X's, 
bodies 1 and 3, and the RCC, body 2. The zone description would be 


and 


A = + 1 + 2 

B = 0R +3-1 0R +3-2 


Notice that the 0R operator refers to all following body numbers until 
the next 0R operator is reached. 
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Fig. 2. Use of 0R operators. 

The geometry must be specified by establishing two tables. 
The first table describes the type and location of the set of bodies used 
in the geometrical description. The second table identifies the physical 
zones in terms of these bodies. The input routine processes these 
tables to put the data in the form required for ray tracing. Because the 
ray tracing routines cannot track across the outermost body, all of the 
zones must be within a surrounding external void so that all escaping 
particles are absorbed. Also no point may be in more than one zone. 

The information required to specify each type of body is as 

follows: 

a. Rectangular Parallelepiped (RPP) 


Specify the minimum and maximum values 
of the x, y, and z coordinates which bound 
the parallelepiped. 











v..v j 




z 



Fig. 3. Rectangular Parallelepiped (RPP). 
b) Sphere (SPH) 

Specify the vertex V at the center and the scalar, 
R, denoting the radius. 



c) Right Circular Cylinder (RCC) 

Specify the vertex V at the center of one base, 
a height vector, H, expressed in terms of its 
x, y, and z components, and a scalar, R, 
denoting the radius. 


^Ymax 
—Zmax 
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Fig. 5. Right Circular Cylinder (RCC). 


d) Right Elliptical Cylinder (REC) 

Specify coordinates of the center of the base 
ellipse, a height vector, and two vectors in 
the plane of the base defining the major and 
minor axes. Presently this body is not 
implemented. 



Fig. 6. Right Elliptical Cylinder (REC). 








w r**r 
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e) Truncated Right Angle Cone (TRC) 

Specify a vertex V at the center of the 
lower base, the height vector, H, expressed 
in terms of its x, y, z components, and two 
scalars, Rj and R2 , denoting the radii of the 
lower and upper bases. 



Fig. 7. Truncated Right Angle Cone (TRC). 
f) Ellipsoid (ELL) 


Specify two vertices, VjandVg, denoting 
the coordinates of the foci and a scalar, R, 
denoting the length of the major axis. 









Specify the vertex V at one of the corners 
by giving its (x, y, z) coordinates. Specify 
a set of three mutually perpendicular 
vectors, ai , with a^ and ^ describing the 
two legs of the right triangle of the wedge. 
That is, the x, y, and z components of the 
height, width, and length vectors are 
given. 



Fig. 9. Right Angle Wedge (WED). 

Box (B0X) 

Specify the vertex V at one of the corners by 
giving its (x, y, z) coordinates. Specify a set 
of three mutually perpendicular vectors, a^ 
representing the height, width, and length of 
the box, respectively. That is, the x, y, 
and z components of the height, width, and 
length vectors are given. 
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Vx, Vy, Vz 

Fig. 10. Box (B0X). 
i) Arbitrary Polyhedron (ABB) 

Assign an index (1 to 8) to each vertex. For 
each vertex, give the x, y, z coordinates. 
Each of the six faces are then described by a 
four-digit number giving the indices of the 
four vertex points in that face. For each face 
these indices must be entered in either clock¬ 
wise or counterclockwise order. 

Y2 
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2.2 


Subroutines 


Descriptions of the combinatorial geometry routines are given 
with logical flow charts depicting their functions. Two M0RSE routines, 
G0MST and L00KZ required significant modifications for use with the 
CG package and descriptions of these routines are also provided. 

2. 2.1 Subroutine G1(S, MA, FPD, L0CREG, NUMB0D, IR0R, 

ml, IR2) 

G1 is the control routine for the combinatorial geometry. On 
one call, it calculates the distance travelled in the present zone, and 
the number IR of the next zone to be entered. Essentially GG is called 
for each body adjacent to the present zone, calculating RIN and R0UT, 
the distances to entry and exit of the body along the trajectory. The next 
zone to be entered is determined by again calling GG to calculate RIN 
and R0UT for each body adjacent to the next possible zone. These next 
possible zones are determined by examining a list cf all the previous 
zones entered on crossing this body. RIN and R0UT are checked against 
the input zone descriptions to determine the correct zone. If it is not 
found in the list of previous zones, all other zones are examined in a 
similar fashion, and when the correct zone is found, it is added to the 
list of previous zones for that body. If the new zone is different from 
the old, G1 returns; otherwise G1 continues tracking until a different 
zone is encountered. One change added to the M0RSE version of G1 is 
that if the distance to the next boundary is greater than the distance to 
scattering, G1 returns without determining the next zone past the bound¬ 
ary, setting the flag MARKG in common 0RGI. 

Called from: G0MST, EUCLID, MESH 

Subroutines called: GG 

Commons required: PAREM, G0ML0C, DBG, 0RGI 
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Variables required: ! 


. i 

XB(3) - starting coordinates of present trajectory, 

WB(3) - direction cosines of trajectory, ^ 

\ ] 

' IR - present zone, 

DIST ’■ - present distance from XB(3), 

DISTO -, distance frbm XB(3) to next scattering point, 

1 ! \ 

NASC - less than zero if this is a! new trajectory, 

'j KL00P - trajectory index, 

PINF - machine infinity, 

MA, FPD, L0CREG, NUMB0D, IR0R, ' 

. IRl, IR2 - locations in blank common used for variable 
dimensioning. 

Variables changed: , 

KL00P - trajectory index incremented if this is a 
new trajectory, 

NASC - next ij>ody intersected by trajectory, 

LSURF - i surface of body NASC crossed at next inter¬ 
section (negative if leaving and positive if 

^ entering NASC), 1 

*DIST - distance from XB(3) to next intersection or 

collision site, . » 

i '• 

MARK - set to 1 if distance to collision! (DISTO) is 

i less than distance to next intersection (other¬ 
wise 0), 1 

S ■ - distance travelled on this call to Gl, 

IRPRIM - zone to be entered on boundary crossing, 

MA'(INEXT) - new zone added to the list of next possible 

zones for body NB0, 

MA(INEX) - location in MA of next item in the next possible 
zone list for body.NB0 (these lists leap-frog 
through the end of the MA array). 
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Significant internal variables: 


NB0 

\ 

R0UT 

RIN 

LRI 
LR0 
! IRP 

Limitations: 

I0UT 


absolute value is body being considered 
white a negative or positive sign indicate^ 
that zone IR or IRP is outside or inside the 
body respectively, 

distance to exit of body NB0 calculated by 
GG, 

distance to entry of body NB0 calculated by 
GG, , 

surface of body NB0 entered by trajectory, 

surface of body NB0 trajectory exits, 

zone being considered as next zone. 


output logical unit is set to 6 in data state¬ 
ment. 
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Subroutine G1 (S, MA, FPD, L0CREG, NUMB0D, IH0R, IR1, IR2) 





























w 















Loop over all remaining regions IRP to determine next zone. 


Loop over all bodies NB0 defining zone IRP. 


Call GG 


NB0 negative?}-►no 


/ Is R0UT<O or\ 

no 

no 

/is RIN< DIST \ 

1 R0UT < DIST 1 



j and j 

\ or RIN >DIST j 


r 

\ DIST < R0UT j 


Found a zone. 
Store it in 
MA(INEXT). 


Error 
Zone not 
defined. 


Dump MA set 
iRPRIM » 0 


(Return 













2.2.2 Subroutine L00KZ (X, Y, Z, MA, FPD, L0CREG, 

NUMB0D, IR0R, NS0RT 

The purpose of this routine is to return the combinatorial geo¬ 
metry zone of point (X, Y, Z) so that tracking can be initialized. The 
coding has been borrowed from the second half of subroutine G1 and 
adapted to determine the zone of a source particle. For efficiency 
L00KZ builds a list of possible source zones to search on future calls. 

If the region is not found on this list, all other zones are examined and 
upon determining the new source zone, it too is added to the list. Notice 
that the starting direction cosines (.8, .6, 0.0) are assumed in L00KZ, 
but may be changed elsewhere. 


Routines called: GG 

Commons required: 0RGI, PAREM, G0ML0C, and DBG. 
Variables required: 


X, Y, Z - coordinates for which a zone is desired, 


MA, FPD, L0CREG, NUMB0D, IR0R, NS0R 

locations in blank common used for variable 
dimensioning. 

Varaibles changed: 

KL00P - trajectory index is incremented, 

NMED - common 0RGI variable set to correct zone 
number, 


NS0R(INEXT) - new zone added to list of possible source 
zones. 


Significant internal variables: 

WB(3) - set to . 8, . 6, 0.0 so that L00KZ need not 

be called with a direction 


* 

CG version. 
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Subroutine L00KZ (X, Y, Z, MA, FPD, L0CREG, NUMB0D, IR0R, 
NS0R) 



RETUR1 







































2.2.3 Subroutine GG (L0CAT, MA, FPD) 


GG is the workhorse of the combinatorial geometry, comput¬ 
ing distances to intersections for all body types. It is called from G1 
or L00KZ to compute distance to entry and distance to exit to a body 
whose location in the MA array is flagged by the argument L0CAT. Each 
time a distance to entry, RIN, and a distance to exit, R0UT, are calcu¬ 
lated for a body they are stored in the MA array together with LRI and 
LR0, the indices of the entry and exit surfaces of the body. Also stored 
at this time is KL00P, the particles trajectory index. On a subsequent 
call to GG for that body, KL00P is checked against the earlier value, 
now L00P. If they are the same, the old values of RIN, R0UT, LRI, 
and LR0 are retrieved so that GG can return immediately. 


If it is necessary to compute a new trajectory a different area 
of coding is entered for each body type to calculate RIN, R0UT, LRI 
and LR0. 


Called from: G1 and L00KZ 

Subroutines called: none 
Commons required: PAREM, GE0M2 
Variables required: 


LOCAT 

KL00P 

PINF 

MA, FPD 


starting locati n in the MA array of integer 
data for the appropriate body, 

trajectory index, 

machine infinity (stored in RIN and R0UT 
when the trajectory misses the body), 

locations in blank common required for 
variable dimensioning. 
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Variables changed: 


RIN 

distance to entrance, 

R0UT 

distance to exit, 

LRI 

- surface of entrance 

LR0 

surface of exit. 

Significant internal variables: 

L 

same as L0CAT, starting location in MA 
array of body data for the appropriate body, 

K 

starting location in FPD array of floating 
point data for the appropriate body. 
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Subroutine GG (L0CAT, MA, FPD) 


^START^ 



Retrieve K, L00P, and 
rrYPE from MA. 






/Has this trajectory previously "N 
V been calculated for this body ? / 


no 


-Yes- 


Set LRI, LR0, RIN 
and R0UT from MA 
and FPD arrays. 






Set RIN or R0UT to DIST when 
one is within roundoff of DIST. 





1r0UT =R0UT*1.OOOO1| 


X If RIN*R0UT ) 


|R0UT = - infinity) 


Write 
"Type not 
implemented' 


Store RIN and R0UT 
in FPD array. 



ETUR1 


RETU) 


















2.2.4 Subroutine GENE (MA, PPD, L0CREG, NUMB0D, IR0R, 

^RIZ, MRCZ, MMIZ, MMCfe, KR1, KR2, il, I0UT, IN) 

Although subroutine J0MEN is used to read in all of the com¬ 
binatorial geometry data except the region volumes, that data is put into 
the proper storage location in blank common by GENI. This is accom¬ 
plished by having J0MIN write the body and zone data on a mass storage 
unit, so that GENE can retrieve that data. GENI assigns the data to 
blank common in the area set aside by J0MEN. For an arbitrary poly¬ 
hedron body, GENI calls ALBERT to handle a portion of the data. In 
addition, GENI outputs the input information concerning the geometry 
data. GENI computes certain geometry dependent data needed during 
the random walk. 


Called from: 

J0MIN 

Subroutines called: 

ALBERT 

Commons required: 

Variables required: 

Blank, G0ML0C, PAREM 

11 

input unit, 

I0UT 

output unit, 

IN 

bulk storage logical unit. 

Variables input: 

Geometry data for zones and bodies is retrieved 
from mass storage unit IN. 

Variables changed: 

MA(I), FPDCE), L0CREG(I), IR0R(I), MRCZ(I), 
MMCZ(I), KR2(I)--see Fig. 13, layout of com¬ 
binatorial geometry storage in blank common. 


i 
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Subroutine GENI (MA, PPD, L0CREG, NUMB0D, IR0R, MRIZ, MRCZ 
MMIZ, MMCZ, KR1, KR2, IN) 


*■ 



RETURN 




















2. 2. 5 Subroutine ALBERT (F, IERR) 

ALBERT is called by subroutine GENI to process the arbi¬ 
trary polyhedron (ARB) body data before storage in the FPD array. The 
ARB body data as read by GENI consists of the coordinates of all corners 
(eight for a six-sided figure), followed by a packed decimal number for 
each side indicating which corner points make up that side. ALBERT 
processes this data, returning a unit normal vector and a minimum 
distance to origin for each plane containing a side of the ARB. An ARB 
can have up to six sides. These unit vectors and distances then replace 
the original ARB data in the FPD array. The number of sides and a 
distance that is characteristic of the ARB's minimum dimension are 
also stored in the FPD array. This minimum distance is later ust... for 
round-off tests. 


Called from: 

GENI 


Subroutines called: 

none 


Commons required: 

Variables required: 

none 


F 

- 

FPD array read from GENI 

Variables changed: 

F 


FPD array changed to contain unit normals 
and distances to origin for each plane. Num 
ber of sides and a minimum distance are 
also stored. 

Significant internal variables: 

X(8,3) 

- 

coordinates of cornei points, 

IS (6,4) 

- 

indicates corners contained in each of the 
possible six surfaces, 

V(3,4) 

- 

edge vectors for a given surface (an inner 


cross product then gives the normal vector). 


29 


JJ 

























JPW5W'> 






i f \ 


2 * 2 * 6 Subroutine PR (Combinatorial Geometr y) 

Subroutine PR is called from various locations in the combi¬ 
natorial geometry package (GENI, Gl, L00KZ and N0RML) whenever 
intermediate or debugging output it required. The amount of geometry 
i data which is out P ut depends on the value of the argument, ‘if this argu¬ 
ment is 1 or 8, all of the geometry data in blank common is printed, 
otherwise only selected variables are output. B^r comparing the argu¬ 
ment value given in the output with the source listing, the geometry data 
at a, given time in the ^execution can be determined. The call to PR is 
initiated by setting the IDBG variable to a nonzero valued Because of 

the large amount of output generated, this option should not be used dur¬ 
ing a normal execution. , ~ 

1 

Called from: GENI, Gl, L00KZ, N0RML 

Subroutine called: none 

1 Commons required: blank, PAREM, G0ML0C, DBG 





Subroutine PR (K) 


( START 

v y 


0: K, XB(I), WB(I), IR 

m, ERPRIM, NASC, LSURF, LRI, LR0, KL00P, 
L00P, ITYPE, N, NUM, L0CAT, ISAVE, 

INEXT, IRP, INEX, LDATA, RIN, R0UT, SMIN, 



K = 1 or 


RETURN 


0: G0ML0C Array 
MA Array 
FPD Array 

L0CREG, NUMB0D, IR0R, MRIZ, MRCZ, MMIZ, 
MMCZ, KRl, KR2, and NS0R Arrays 
Region Volume Array 















2.2.7 Subroutine N0RML (MA, FPD, L0CREG, NUMB0D )* 

The purpose of subroutine N0RML is to return a unit vector 
to a combinatorial geometry body NASC at point XB + D 1 * WB which 
must be on the surface LSURF of body NASC. This unit vector is useful 
either for albedo scattering or boundary crossing flux estimates. The 
sign of the unit vector is chosen so that wB • n is negative meaning that 
the unit vector will point against the particle direction. 

Called from: ALBD0 

Commons required: N0RMAL, G0ML0C, PAREM 

Variables required: 

NASC - body number particle is on, 

LSURF - surface of body NASC r 

XB(3) ) 

WB(3) > - trajectory parameters, 

DIST ) 

MA, FPD, L0CREG, NUMB0D 

locations in blank common used for variable 
dimensioning. 

Variables changed: 

UN(3) - direction cosines of unit vector stored in 

common N0RMAL. 

Significant internal variables: 

XP(3) - set to XB + DIST * ifiTJ, 

X(3) - usually used to relate the intersection point 

to a body centered coordinate system, 

H(3) - usually used as a body orientation, 

L - body location in MA array, 

K - body location in FPD array. 

It 

CG version. 
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2.2.8 Subroutine r~VLIN (FPD, MRIZ, VN0R, IV0PT, NIR, 

niz, ii, 10; ----— 

This routine is called by J0MIN to read in or to calculate the 
volume of each region in the geometry. The four options available are 
(1) to set each volume to 1., (2) to calculate the volumes for concentric 
spheres, (3) to calculate the volumes for slabs (not coded at present), 
or (4) to read in the volumes for each region from cards. The volumes 
are stored in blank common and are only used by the track length and 
collision density estimators. 


Called from: 

J0MIN 


Subroutines called: 

none 


Commons required: 

Variables required: 

none 


FPD 

- 

array containing zone description data, 

MRIZ 

- 

array which relates M0RSE region to input 
zones, 

VN0R 

- 

array containing the M0RSE volumes, 

IV0PT 

- 

options for determining volumes. 

NIR 

- 

number of regions, 

NIZ 

- 

number of zones, 

11 , 10 

- 

input and output logical units. 


Variables changed: 

Variables in blank common starting at KV0L. 




















2.2.9 Subroutine J0MIN (NADD, II, IQ) * 

The input of the geometry data is controlled by the J0MIN 
subroutine, which performs the following tasks: 

• Reads all geometry input data except the region volumes. 

• Writes the body and zone data on a mass storage unit 
(I0UT-16). 

• Determines the length of all geometry arrays. 

• Calculates the beginning location in blank common of 
geometry arrays. 

• Initializes geometry arrays. 

• Calls the GTVLIN subroutine which returns region 
volumes. 

Since combinatorial geometry input data is dynamically allocated to con¬ 
serve storage area, it is stored temporarily on a mass storage device. 
This allows the core storage requirements to be determined. Hence, 
much of the coding in J0MIN is similar to GENI, which reads the data 
on the mass storage device and puts it into blank common. 

Called from: INPlfol 

Subroutines called: GENI, GTVLIN 
Commons required: B^ank, G0ML0C, PAREM 
Variables required: ^ 

4 

All variable^' in G0ML0C, II, 10 

4 

Variables input: : 

IV0PT, IDBp, NAZ, MRIZ(I), MMIZ(I) 

A 

Variables changed: A 

A 

All variable^ in G0ML0C, IV0PT, IDBG, MRIZ(I), MMIZ(I). 

A 


CG version. 








Important internal variables: 


I0UT 

NAZT 


mass storage unit, 

total number of adjacent zones summed 
all zones. 


over 






Subroutine J0MIN (NADD, 11, 10) 



: Cards CGD and CGE. 


Call GENI 


Calculate number of regions (NIR) and total 
storage required for geometry data (NADD) 


Zero the NS0R array. 


Call GTVLIN. 


















2 . 2. 10 Subroutine GTV0L (MXREG, GN0R) 

This routine is called by ENDRUN to calculate the reciprocal 
of the volume of each region. It is geometry independent, since it uses 
information determined in GTVLIN. 

Called from: ENDRUN 

Subroutines called: none 
Commons required: G0ML0C, blank 
Variables required: 

KV0L - index from common G0ML0C, 

MXREG - number of regions in the geometry. 

Variables changed: 

GN0R - reciprocal of the volume ior p ion. 
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2.2.11 Subroutine G0MST (TSIG, MARK) * 

Any boundary crossing between the present and next collision 
sites are determined by calling the combinatorial geometry routine Gl. 
Before the Gl call, combinatorial geometry variables in common PAREM 
are initialized, and after the call NUTR0N variables are updated. 

WATERS is summed in blank common as a track length flux estimator. 
MARK is set to -1 for an external void and -2 for an internal void. 


Called from: 

NXTC0L 

Subroutines called: 

Gl 


Commons required: 

NUTR0N, AP0LL0, PAREM, 0RGI, G0ML0C 

Variables required: 



X, Y, Z, U, V, W, NMED see common NUTR0N for definitions. 

BLZNT 

- 

value of IR from last track or from L00KZ, 

XB(3) 

- 

starting coordinates of present trajectory, 

WB(3) 

- 

trajectory direction cosines, 

MARK 

- 

flag to indicate type of trajectory. 

ETATH 

- 

distance to be travelled in cm if the flight 
remains in the same media, 

DIST 

- 

present distance from XB(3). 

Variables changed: 



X, Y, Z 

- 

endpoints of flight, 

DISTO 

- 

distance from XB(3) to next collision site, 

ETAUSD 

- 

mean free paths travelled on this call to 
G0MST, 

ETATH 

- 

cm travelled on this call to G0MST, 


* 

CG version. 
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MARK - flag indicating type of termination of flight, 

0 - normal boundary crossing, 

1 - flight within one medium, 

-1 - particle escaped, 

-2 - particle entered an interior void, 

NMED and 

NREG - medium and region of end point, 

combinatorial geometry region of end point. 


BLZNT 






















2.2.12 Subroutine EUCLID (MRK, XI, Yl, Zl, X2, Y2, Z2, P1P2, 
IG ) ARG, NT, MEDIUMTIBLZ, NREGN )* 

Subroutine EUCLID is used to determine the number of mean 
free paths between two points. However, if the argument NT is non¬ 
zero only the mean free 'paths to the first media boundary will be calcu¬ 
lated. The routine works much like NXTC0L and G0MST. After ini- 

tializihg the combinatorial geometry variables, the particle is tracked 

' \ \ 

from point one to point two. Total cross sections are determined by 
calling NSIGTA and distances to next interfaces are calculated by call¬ 
ing Gl. On an initial call to EUCLID, MRK should be 1. If NT is non¬ 
zero, MRK should be 0 on successive calls for the same trajectory, so 
that new trajectory parameters are not initialized. , This will be handled 
automatically if the calling routine does not change MRK. 


Called from: I 

GETETA, RELC0L 

Subroutines called: 

Gl, 

NSIGTA 

Commons required: 

Variables required: 

PAREM, 0RGI 

\ 

1 

MRK 

- 

set to 1 for initial call of a trajectory; 

XI, Yl, Zl 

- 

coordinates of starting point, 

X2, Y2, Z2 

- 

coordinates of end point, 

' P1P2 

- 

distance between starting and end points. 

IG ' 

- 

energy group index, 

NT 

- 

0 for total mean free paths, 4 0 for mean 
free paths between intersection points, 

MEDIUM 

- 

NMED of point one, '« 

NREGN 

- 

NREG of point one, i 

IBLZN 

- 

IR of point one. 


CG version. 
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Variables changed: 
MRK 


1 for a flight reaching the end point, 

0 for a flight crossing a medium boundary 
(NT 4 0 only), 

-1 for a flight escaping the system, 

-2 for a flight encountering an internal void 
(NT 4 0 only), 

Xl, Yl. Z1 - returns boundary intersection point if 

NT 4 0, 

ARG - negative of number of mean free paths, 

NT - if NT 4 0 on input, will return as -1 if an 

escape occurs, 

MEDIUM - medium number of end point. 

Significant internal variables: 

ETA ~ distance remaining to point two in cm, 

ETAUSD - distance travelled on last call of Gl. 





















2. 2.13 Subroutine ENDRUN 


This subroutine is called at the end of a run. The weight of 
the particles at collision and the track length of the particles has been 
summed over all collisions within the region for all particles in the run. 
Thus, a normalization by the number of particles in the run is required. 

In addition, the volume of the region and the total cross section for each 
region and energy group is used to determine the average fluence over 
the regions for the collision density estimate. For the track length esti¬ 
mate the same coding is used with the total cross section set to 1. Note 
that a region cannot contain more than one medium. Since storage al¬ 
locations are flexibly dimensioned, NEX in the SAMB0 input must provide 
for MXREG + 2 arrays in Blank Common to be used by ENDRUN. A 
check is made and if this condition is not met, no output from ENDRUN 
is obtained. 

There are many manipulations in ENDRUN with the energy 
indexes. These manipulations are required in order to treat the six 
major options in M0RSE. Appendix A gives an example of the six prob¬ 
lems. 

Called from: NRUN 

Subroutines called: GTV0L, NSIGTA 
Commons required: blank, PDET, USER 
Variables required: 

NMED (REGI0N) - must be set in data statement, 

GN0R - reciprocal of the volume for each 

region. 

Many variables in blank common. 
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Significant internal variables: 

FNT - total number of source particles in the run, 

MXREG - the number of geometry regions, 

LIM1 - the group number for the last primary 

particle group, 

EB0T1 - lowest primary particle group energy, 

EB0T2 - lowest secondary particle group energy, 

NGP3 - the group number for last secondary group, 

NQP1 - the group number for last primary group 

(negative if gamma-ray only problem), 

IFLG - flag to indicate which fluence estimate is 

currently being output. 0-collision density, 
1-track length. 







Subroutine ENDRUN 



(Collision density?) -Yes—► f~Set titles 1 

No IFLG=0 


1 Set titles.) 


|Calculate lowest neutron and gamma-ray energy, j 


i___ 

(Calculate group widths, AE.J 


' 

(Adjoint problem?)-Yes-^ Set AE = 1| 


IADJ > 0 

No 


(Call OTV0L calculates 1/volume of each region. | 


0 







m 


I < NMGP 


Arc there more primary 
particle groups? 


to- (Are there secondary particles? 


Yes 

NMTG / NMGP 


Call NSIGTA to find total cross section 


Store fluence estimate in blank common, 
sum over energy (secondary). 
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(Adjoint problem ?)—Yes—*-| 

0: Heading 


IADJ>0 I 

for table. 


No 



|0: Heading for table.) 




| 0: Results | 





[Calculate Ee>(E) for each primary group.) 




- {Are there secondary particles?) 


Yes 

NGPQTN+NGPQTG > C 


Calculate Ecn(E) for each secondary group. 
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III. MODIFICATIONS FOR UNIVAC-1108 

In converting the M0RSE-SAMB0-M0RSEC code frorr the 
IBM 360 to the UNIVAC-1108, there were a considerable number of 
changes that had to be made. These changes were connected mostly 
with the local library functions of the IBM-360 and an attempt was made 
to not add local UNIVAC-1108 library functions. The only reduction in 
the code's capability is in the diagnostic module where the ability to 
look at a bit string and decide whether it was a floating point or integer 
number, convert it to hollerith and output it, does not exist in this version. 
However, the diagnostic module is mostly machine independent and still 
has the capability of writing out values from common and parts of blank 
common. Subroutines BNKHLP and HELPER are dummy routines. A 
diagnostic routine for the combinatorial geometry package has been added. 

Routines which are still machine dependent (involving coding 
other than word lengths, data statements and equivalences) are DATE 
and TIMER. Subroutine DATE calls ERTRAN and uses DEC0DE. Sub¬ 
routine TIMER c. !ls CPUTIM that assumes that the time is returned in 
units of microseconds. 

Subroutine READSG is still machine dependent in its need to 
read each card twice and subroutine ERR0R has been added—its only 
function is to call EXIT. Subroutine INPUT now calls INPUT1 and 
INPUT2 with this change, permitting a better use of overlay. The 
first six variables in NXTRA from common AP0LL0 are used. 

The layout of blank common is given in Figs. 12-18 for the 
different modules of M0RSE. Definitions of variables in blank common 
are given in Table HI which defines the indexing scheme. Variables in 
the labelled commons are defined in Tables IV - XVI. 

Preceding page blank 
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A description of the energy indexing scheme in M0RSE- 
SAMB0-M0RSEC is given in Appendix A. 

3.1 Built-In Fluence Estimators 

There are two fluence estimators that are "built-in” the 

UNIVAC 1108 combinatorial geometry version of M0RSE. The addition 

of the collision density estimator to M0RSE has been discussed pre- 
(7) 

viously v and the track length per unit volume estimator has been added 
to take advantage of a feature of the combinatorial geometry. The re¬ 
sulting fluence is averaged over geometry regions and if no regions are 
used only the spectrum averaged over the whole system is available. 


3.2 Collision Density Fluence Estimator 

An estimate of the fluence averaged over a region specified 
by the geometry input may be obtained from the weight that has been 
summed over all collisions in the region by 


VV - 


V TVi 


where 


V E i> 


= fluence in region V at energy E., 

- number of particles of energy E. scatter¬ 
ing in region V, 

= weight of the n n particle of energy E. 


NITS*NSTRT = 


S TV (E i> 


scattering in region V, 
total number of particles, 
volume i/ region V, 
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tota 1 cross section for groups E. in region 
V. 1 









3.3 


Track Length per Unit Volume Estimator 

Another estimate of the fluence averaged over a region may 
be obtained from the track length of the neutron or photon in that region. 
The track length calculated by the combinatorial geometry package is 
stored in a manner similar to the collision weights. The fluence esti¬ 
mate may be obtained by 



N Vi = number of trajectories of energy E. in region V, 


1L 

w n = weight of particle on n trajectory, and 

tfl 

L = track length of n trajectory in region V at energy 

n E. . 

1 


Subroutine ENDRUN is called to process the weights and 
track lengths to obtain estimates of the fluence. In an adjoint problem, 

the importance is calculated. The fluence estimate has units of particles/ 

2 2 
eV/cm /source and in an adjoint problem the units are importance/cm / 

source. The total fluence is integrated over energy from the primary 
and secondary energy groups to obtain an integral value for each region. 
The only input required for this routine is a data statement giving the 
medium number for each region; there is a limitation of only one medium 
^gion. In addition, one of the several options for inputing or calcu¬ 
late .he volume of each region must be utilized. 

The major disadvantage with using these estimators is that 
no estimate of statistics is available. The use of this estimator is in¬ 
expensive, however (output routine is called only once), and is meant 
to be used as a supplemental estimator. 
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Starting 

Location 

Mnemonic Variable 
Name 

Length 


KNER 



VEL 

! 4*NMTG 


FS 

L0CWTS 

nps 


) 


Current Wei 
Standards 

>ht 

| 3+MGPREG 


Path 

| MGPRFG 


SPLIT and R.R. 
Counters 

J 8*MGPREG 

L0CFWL_.. 

Initial Weight 
Standards 

J 3*MGPREG 


Current FWL0 

J 2*MXREG 


initial FWL0 

L.0CEPR 

GW 1.0 


[ NMTG*MXREG 

LOCNSC 

EPRB 

} NMTG*MXREG 


Scattering; 

Counters 

} 10+NMTG+MXREG 

L0CrSN_.. 

NSCA 

| MEDIA 


FISH 

| NMTG*MEDIA 


FSE 

j NMTG* MEDIA 

NGE0M 

GMGN 

| NMTG+MEDIA 

NGLAST_* 

Geometry Data 

j See Figs. 13-15 for 
details. 

NS1GL 

Permanent Cross 
Sections 

| See Fig. 16 for details. 

NM)>SEC . _ 

Temporary 

Cross 

Sections 

Particle 

Bank 

1 14*NMOST see Fig. 17 
j for details. 

NLAST 

Fission 

Bank 

\ 7 "N MO ST IF MFISTP>0 
j see Fig. 18 for details. 


User Area 

1 NLEFT see Fig. 19 for 
j details. 


Fig. 12. General layout of blank common. 
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I 


Starting 

Location 


information 


Size 


NGE0M=NADD 

KMA 

KFPD 

KLCR 

KNBD 

K30R 

KRIZ 


Length Df geometry array 

MA ' 

Integer array j 

FPD 

Floating point array 
L0CREG 

Indices to correlate MA array 
data wi n code zone data ■ 

NU1^B0D 

Number of bodies for each 
code zone 

IR0R 

Indices to correlate input zone 
to code zone : 

- [ - 1 - 

MRIZ 

Indices to correlate M0RSE 
region to input zone 


KRCZ — 


KMIZ 

KMCZ 

KKR1 

KKR2 j 

KNSR 

KV0L 

NGIAST 


MRCZ 

Indices to correlate M0RSE 
region to code zone 

NMIZ 

Indices to correlate M0RSE 
media to input zone 


NMCZ 

Indices to correlate M0RSE 
media to code zone 

KRlj ’ 

Indices to correlate first code 
zone to ‘ nput zones 

KR2 

Indices to correlate last code 
zone to input zone 


NS0R 

Indi-.cs of code zones in which 
source particles have been found 

VN0R 

i Volume of each M0RSF region 


LTMA 


LFPD 


MUMR 

I 


NUMR 


NUMR 


IRTRV 


NUMR 


IRTRV 


NUMR 


IRTRV 


IRTRV 


NUMR 


NIR 


Fig. 13. Layout of combinatorial geometry data in bla *xiv common. 
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Figure 15. 








Position in 
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Location _ 

1RSG = NLAST 

= NGLAST+1 

INGS 

INSG 

IDEL 


_ information 

List of Mixing Table 


Index to Egg (Primary) 

Index to Egg (Secondary) 

List of Element I.D. Numbers 


INNN 

IGGG 

ISTART 

JSCC0G 

INAB0G 

IGAB0G 

IFP0RG 

INUS 

INGN 

INGNP 

IFNGT> 

IFSP0G 

IDSG0G 

IPRBNG 

IPRBGG 

ISCANG 

ISCAGG 

ISP0RG + ISTART 
If ISTAT > 0 ISP0RT 


NLAST or NNIC 
NNIC 


Number of Downscatters for each 
primary group 

Number of Downscatters for each 
secondary group _ 


s 


V e t 


v£ f /S T 


l . „ 


Ej^/E^ not used at present 
E not used at present 


'N-r 
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g -B 
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Repeat for next medium 


l'i Coefficient Primary 
Pj Coefficient Secondary 
Repeat for P^ Coefficient 

R epeatTor next medium 
Index to point cross sections for 
each multigroup boundary 


Size _ 
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NGG 

NELEM*NC0EF 
if IXTAPE > 0 
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NGG 
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NDSNGG* NSCT 

NDSNGP* NSCT 
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NDSNGG 

INFP0G*(NC0EF-2) 

INF P0G*(NC0EF-1) 
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I06RT>O; otherwise 
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Location 


Information 


Size 


NXSECT(l) 


NXSECT(J) 


NLAST 


Point cross 
sections for 
Medium 1 

total 

scattering 

NPT(1)*(NEGPS-1)* 

NXPM 

Point cross 
sections for 

total 

NPT(J)*(NEGPS-1)* 

Medium J 

scattering 



NXPM 


Fig. 16. Layout of permanent cross sections in blank common. 
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Location 


Variable 


';N0 = NSIGL 


NN0 + 14 


IG 

U 

V 

w 

X 


Y 

z 

WATE 

AGE 

BLZNT 

NAME 

NAMEX 

NMED 

NREG 


Variables are from 
NUTR0N common, 
see Table VI for 
definitions. 


Repeat for Particle 2 


NN0 + 14*NM0ST = NLAST 

= NFISBN IF MFISTP > 0 


Fig. 17. Layout of particle bank in blank common. 








Location 

Variable 

NFISBN 

X 


Y 


Z 


WATEF 


AGE 


IG 


NAMEX 

NFISBN + 7 

Repeat for particle 2 

NFISBN + 7*NM0ST = NLAST 

• 


Fig. 18. Layout of fission bank in blank common. 




Location Labels 


Mnemonic 
Variable Name 


Length 





NLAST*+1' .. . .~ 

L0CRSP 

LABELS 

20*NRESP 



RESP 

NRESP*NMTG 

L0CXD 


EXTR 

NEX*NMTG 



XD 

YD 

ZD 

RAD 

T0 

FACT 

6*ND 

L0CIB 


EXTR 

NEXND*ND 

L0CC0 

IB 

EP 

DELE 

3*NE 

L0CT 

C0S 

NA 

L0CUD 

T 

DELT 

2*ND*NT 

L0CSD 

UD 

SUD 

SUD2 

3*ND*NRESP 



SD 

SSD 

SSD2 

3*ND*NRESP 

L0CQE 


SUD & SSD Units 

20 





NOTE: NLAST must not be changed after SC0RIN is called. 
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Location Labels 

Mnemonic. 
Variable Name 

Length 


QE 



SQE 

3*NE*ND 


SQE2 


L0CQT__ 

SQE Units 

20 


QT 



SQT 

3 *NT*ND*NRESP 


SQT2 


L0CQTE 

SQT Units 

20 


QTE 

SQTE 

SQTE2 

3*NT*NE*ND 

L0CQAE_J 

SQTE Units 

20 


QAE 

SQAE 

SQAE2 

3*NA*NE*ND 

LMAX_, 

SQAE Units 

20 





Fig. 19. Layout of user area (analysis) in blank common. 
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Indexing of Random Walk Blank Common Arrays 
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Response function. 





TABLE IV. (Cont'd.) 








TABLE V. 


Definition of Variables in AP0LL0 Common 

Variable Definition 

AGSTRT Input starting age of source particle. 

DDF Starting particle weight as determined in S0RIN. 

DEADWT(5) The summed weights of the particles at death. The 

four deaths are: Russian roulette, escape, energy, 
and age limit. DEADWT(5) is unused. 

ETA Mean-free-path between collisions. 

ET/-.TH Distance in cm to the next collision if the particle 

does not encounter a change in total cross section. 

ETAUSD Flight path in m. f. p. that has been used since the 

last event. 

UINP, VINP, WINP Input direction cosines for source particle. 

WTSTRT Input starting weight. 

XSTRT, YSTRT, 

ZSTRT Input starting coordinates for source particle. 

TCUT Age limit at which particles are retired. 

XTRA(IO) Not used. 

10, II Output and input logical units. 

MEDIA Number of media for which there are cross sections. 

IADJM Switch indicating an adjoint problem if > 0. 

ISBIAS Switch indicating that source energy distribution is 

to be biased if > 0. 

IS0UR Input source energy group if > 0; otherwise, S0RIN 

is called to read input spectrum. 

ITERS Number of batches still to be processed in the run. 

ITIME Not. used. 

ITSTR Switch indicating that secondary fissions are to be 

the source for the next batch if > 0. 
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TABLE V. (Cont'd.) 


Variable 


Definition 


L0CWTS 

L0CFWL 

L0CEPR 

L0CNSC 

L0CFSN 

MAXGP 


M/L^IM 


MEDALB 

MGPREG 

MXREG 

NALB 

NDEAD(5) 

NEWNM 

NGE0M 

NGPQT1* 

NGPQT2* 

NGPQT3* 


Starting location in blank common of the weight standards 
and other arrays MGPREG long. 

Starting location in blank common of the fission weights. 

Starting location in blank common of the energy-biasing 
parameters. 

Starting location in blank common of the scattering 
counters. 

Starting location in blank common of the fission and 
gamma-generation probabilities for ea-.h medium and 
group. 

Maximum number of energy groups for which there are 
weight standards or path-length stretching parameters. 

The elapsed clock time at which the problem is termin¬ 
ated. 

Medium number for the albedo medium. 

Product of number of weight standard groups (MAXGP) 
and regions (MXREG). 

Maximum number of regions in the system. 

An index indicating that an albedo scattering has occurred 
if > 0. 

Number of deaths of each type (see DEADWT). 

Name of the last particle in the bank. 

Location of first cell of geometry data storage in blank 
commc 

The lowest energy group (largest group number) for 
which primary particles are to be followed. 

The number of primary particle groups. 

The lowest energy group (largest group number) for 
which any particle is to be followed. 


* 

See Appendix A for diagram of energy group structure. 
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TABLE V. (Cont'd.) 


Variable 


Definition 


NGPQTG 

Number of energy groups of secondary particles to be 
followed. 

NGPQTN* 

Number of energy groups of primary particles to be 
followed. 

NITS 

Number of batches per run. 

NKCALC 

The first batch to be used for a k calculation. If 0, k 
is not calculated. 

NK3LL 

An index to indicate that Russian roulette is to be 
played if > 0. 

NLAST 

The last cell in blank common that was used by the 
cross-section storage or is set aside for banking. 

NMEM 

The location of the next particle in the bank to be pro - 
cessed. 

NMGP* 

The number of primary particle groups for which there 
are cross sections. 

NM0ST 

The maximum number of particles that the bank can hold. 

NMTG* 

The total number of energy groups (both primary and 
secondary) for which there are cross sections. 

N0LEAK 

An index which indicates that nonleakage path-length 
selection is to be used if >0. 

N0RMF 

An index to indicate that the fission parameters are to 
be renormalized if >0. 

NPAST 

An index to indicate that the exponential transform is 
to be used if > 0. 

NPSCL(13) 

An array of counters of events for each batch: 

1. sources generated 

2. splittings occurring 

3. fissions occurring 

4. gamma rays generated 


* 

See Appendix A for diagram of energy group structure. 
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TABLE V. (Cont’d.) 

Variable Definition 


NQUIT 

NSIGL 

NS0UR 

NSPLT 

NXTRA(l-5) 

NXTRA(6) 

NXTRA(7-10) 


5. real collisions 

6. albedo scatterings 

7. boundary crossings 

8. escapes 

9. energy cutoffs 

10. time cutoffs 

11. Russian roulette kills 

12. Russian roulette survivors 

13. gamma rays not generated because bank was full. 
Number of runs still to be processed. 

Starting location of the bank in blank common. 

An index input to indicate that fissions are to be the 
source for future batches. 

An index to indicate that splitting is to be considered 
if >0. 

Used to separate INPUT1 and INPUT2. 

NGPREG calculated in M0RSE and used in G0MST as 
an index for track length estimation. 

Not used. 




TABLE VI. 


Definition of Variables in NUTR0N Common 


Variable 

Definition 

NAME 

Particle's first name. 

NAMEX 

Particle's family name. (Note that particles do not 
marry.) 

IG 

Current energy group index. 

IG0 

Previous energy group index. 

NMED 

Medium number at current location. 

MED0LD 

Medium number at previous location. 

NREG 

Region number at current location. 

U, V, W 

Current direction cosines. 

U0LD, V0LD, 
W0LD 

Previous direction cosines. 

X, Y, Z 

Current location. 

WATE 

Current weight. 

0LDWT 

Weight at previous collision. 

WTBC 

Weight just before curreni collision. 

BLZNT 

Current block and zone number (packed) for 05R geo¬ 
metry, and current zone number IR for combinatorial 
geometry. 

BLZ0N 

Previuus block and zone number (packed) for 05R geo¬ 
metry, and previous zone number IR for combinatorial 
geometry. 

AGE 

Current age. 

0LDAGE 

Previous age. 






TABLE Vn. 


Definitions of Variables in USER Common 


Variable Definition 


AGSTRT 

WTSTRT 

XSTRT 

YSTRT 

ZSTRT 

DFF 

EB0TN 

EB0TG 

TCUT 

10 

II 

IADJM 

NGPQT1 

NGPQT 

NGPQT3 

NGPQTG 

NGPQTN 

NITS 

NLAST 

NLEFT 

NMGP 

NMTG 

NSTRT 


Initial chronological age to be assigned to source 
particles. 

Initial weight to be assigned to source particles. 

Initial x position to be assigned to source particles. 

Initial y position to be assigned to source particles. 

Initial z position to be assigned to source particles. 

Normalization for adjoint problems. 

Lower energy boundary (eV) of last neutron group 
(group NMGP). 

Lower energy boundary (eV) of last gamma-ray group 
(group NMTG). 

Chronological age limit. 

Logical unit for output. 

Logical unit fcr input. 

Adjoint switch (>0 for adjoint problem). 

Problem-dependent energy group limits (see Appendix 
A). 

Group number of lowest energy gamma-ray group to be 
treated. 

Group number of lowest energy neutron group to be 
treated. 

Number of batches to be run. 

Last cell in blank common used by random walk package. 
Number of cells in blank common available to user. 
Number of primary particle energy groups. 

Total number of energy groups. 

Number of source particles for each batch. 







TABLE VIII. 

Definitions of Variables in Common G0ML0C 


Variable 


Definition 


KMA 

Starting location for the array MA containing integer 
data for each code zone. 

KFPD 

Starting location for the array FPD containing real data 
for each code zone. 

KLCR 

Starting location for the array L0CREG(I) which contains 
the starting location in the MA array for the I th code 
zone data. 

KNBD 

Starting location for the array NUMB0D(I) which con¬ 
tains the number of bodies for the I*h code zone. 

KI0R 

Starting location for the array IR0R(I) which contains 
the index of the corresponding input zone for- the I th 
code zone. 

KKIZ 

Starting location for the array MIRZ(I) which contains 
the index of the M0RSE region corresponding to the I th 
geometry input zone. 

KRCZ 

Starting location for the array MRCZ(I) which contains 
the index of the M0RSE region corresponding to the I^ 1 
geometry code zone. 

KMIZ 

Starting location for the array MMIZ(I) which contains 
the index of the M0RSE media corresponding to the I** 1 
geometry input zone. 

KMCZ 

Starting location for the array MMCZ(I) which contains 
the index of the M0RSE media corresponding to the I** 1 
code zone. 

KKR1 

Starting location for the array KR1(L) contains the first 
code zone which was made from the input zone. 

KKR2 

Starting location for the array KR2(L) contains the last 
code zone which was made from the input zone. 

KNSR 

Starting location for array NS0R which contains the code 
zones in which source particles have been found. 
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TABLE Vm. (Cont’d.) 


Variable 

Definition 

KV0L 

Starting location for the array VN0R(I) which contains 
the volume for M0RSE region (I). 

NADD 

Starting location for the geometry data length and 
changed in J0MIN to the total number of words re¬ 
quired for geometry data. 

LDATA 

Length of the integer data in the MA array excluding 
the words set aside for zone search information. 

LTMA 

Total length of the MA array. 

LFPD 

Length of the FPD array 

NUMR 

Number of code produced zones. 

IRTRJ 

Number of input zones. 

NUMB 

Number of bodies. 

NIR 

Number of M0RSE geometry regions. 









TABLE IX. 


Definitions of Variables in Common L0CSIG 


Variable 

Definition 

ISTART 

Starting location for the total cross-section vector for 
the first medium. 

ISCC0G 

Starting location i„ . the scattering cross-section veccor 
for the first medium. 

INAB0G 

Starting location for the non-absorption vector for the 
first medium. 

IGAB0G 

Starting location for the gamma-ray production vector 
for the first medium. 

IFP0RG 

Starting location for v S f , the fission neutron production 
vector for the first medium. 

IFNGP 

Starting location for the primary-secondary transfer 
probability matrix. 

IFSP0G 

Starting location of the primary downscatter probability 
matrix. 

IDSG0G 

Starting location of the secondary downscatter prob¬ 
ability matrix. 

IPRBNG 

Starting location of the primary scattering angle prob¬ 
ability matrix. 

IPRBGG 

Starting location of the secondary scattering angle prob¬ 
ability matrix. 

ISCANG 

Starting location of the primary scattering angle matrix. 

ISCAGG 

Starting location of the secondary scattering angle 
matrix. 

ISP0RG 

Size of storage needed for each medium, not including 
Legendre coefficients. 
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TABLE IX. (Cont'd.) 

Variable Definition 


ISP0RT* 

INPBUF 

ISIG0G 

INFP0G* 

IABS0G 

IT0TSG* 

NGP 

NDS 

NGG 

NDSG 

INGP 

INDS 

NMED 

NELEM 


Starting location for temporary storage of downscatter 
matrix. 

Starting location for temporary storage of the P^ table. 

Starting location for temporary storage of total cross 
section for element 1. 

Starting location for temporary storage of vLf for 
element 1. 

Starting location for temporary storage of downscatter 
matrix for P^ coefficients (primary groups, element 1). 

Total storage required for temporary storage. 

The number of primary groups to be treated. 

Number of downscatters for NGP (usually equal to NCP). 

Number of secondary groups to be treated. 

Number of downscatters for NGG (usually equal to NGG). 

Number of groups for which cross sections are to be 
inp it. 

Number of downscatters for the INGP groups. 

Number of media for which cross sections are to be 
stored--should be same as MEDIA as read on Card B 
of M0RSE input. 

Number of elements for which cross sections are to be 
read. 


If Legendre coefficients are to be restored, then: 

INFP0G - Redefined by JNPUT as number of locations required 
for each coefficient (both primary and secondary). 

IT0TSG - Redefined by JNPUT as total storage required for all 
coefficients for each medium. 

ISP0RT - Redefined by JNPUT as starting location of Pj co¬ 
efficient for primary groups for medium 1. 
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TABLE IX. (Cont'd.) 


Variable Definition 

NMIX Number of elements times density operations to be per¬ 

formed. 

NC0EF Number of coefficients, including Pq . 

NSCT Number of discrete angles (usually NC0EF/2j nte g ra j). 

NTS Number of downscatters for combined primary and 

secondary groups (usually equal to NTG). 

NTG Total number of groups (primary + secondary) = NGP 

f NGG. 

NDSNGP The number of locations needed for tne downscatter 

matrix for the primary particles. 

NDSNGG The number of locstions needed for the downscatter 

matrix for the seconc \ry particles. 

IADJ Same as IADJM. 

NME Indicator for stripping gamma rays from a coupled 

neutron gamma-ray cross-sectior ;.et--set equal to 
number of neutrons groups + 1. 

L0C Same as L0CEPR. 

INGS Starting location of the indices f< r starting location of 

the downscatter vector for each group for primary 
particles. 

INSG Starting locat’- of the indices for starting location of 

tne downscatter vector for each group for secondary 
particles. 

II, 10 Input and output logical unit numbers. 

KKK A running index of the number of cross sections that 

have already been read in (used in checking the element 
numbers obtained from tape). 

IXTAPE Logical tape number of the multigroup cross-section 

tape if > 0, or logical tape number of the processed 
cross-section tape if < 0. 

IDEL Starting location for element identifiers which determine 

the element cross sections to be read from tape. 


86 









TABLE IX. (Cont'd.) 


Variable 

Definition 

ITEML 

Amount of storage for primary and secondary group 
downscatters per element. 

ITEMG 

Starting location for temporary storage of downscatter 
matrix for Pl coefficients (secondary groups) for 
element 1. 

IRSG 

Starting location of the mixing parameters. 

IRDSG 

Switch to print the cross sections and to test the card 
sequence as they are read if > 0 (test card sequence 
only if = 0, and does neither if <0). 

ISTR 

Switch to print cross sections as they are stored if 
>0. 

IPRIN 

Switch to print angles and probabilities if >0. 

IFMU 

Switch to print intermediate results of n's calculation 
if > 0. 

IM0M 

Switch to print moments of angular distribution if > 0. 

IDTF 

Switch to signal that input format is DTF-IV format if 
>0; otherwise, ANISN format is assumed. 

ISTAT 

Flag to restore Legendre coefficients for next-flight 
estimates if >0. 

IPUN 

Switch to print results of bad Legendre coefficients if 
>0. 

NUS 

Number of groups of upscatter. 

NGN 

Not use a presently. 

IHT 

Location of S-p in the cross-section table, currently 
set to 3 in XSEC. 

INUS 

Starting location for the upscattering vector for the 
first medium. 

INUSN 

Number of impossible first moments found in the cross 
sections. 

INGN 

Starting location for the photoneutron production vector 
for the first medium (not used at present--same as 
INUS). 
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TABLE IX. (Cont'd.) 

Definition 

Starting location for the secondary-primary transfer 
probability matrix (not used at present--same as INUS) 

Starting location for the array of the number of down- 
scatter groups for each primary group. 

Starting location for the array of the number of down- 
scatter groups for each secondary group. 


Variable 

INGNP 

INNN 

IGGG 
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TABLE X. 

Definitions of Variables in Common GTSC1 


Variable 

Definition 

NPT(16) 

Number of points per supergroup for each medium. 

(Only 16 media are allowed.) Redefined at end of 
GTSCT as the total storage required for each cross 
section for each medium. 

NXSECT(17) 

Starting location of point cross sections for each 
medium. 

NDSGP 

Total storage required per supergroup for all media. 

11, 10 

Input and output logical unit numbers. 

I06RT 

Logical tape number of 06R point cross-section tape. 

IGQPT 

Last M0RSE multigroup for which point cross sections 
will be used NMGF> 

NEGPS 

Number of supergroup boundaries. 

ESPD(IOO) 

Energy boundaries (in eV) of the supergroups (only 100 
supergroups are allowed). 

NNIC 

Starting location of the point cross sections. 

NXPM 

Number of cross sections per medium 

= 1 if total only 

= 2 if total + scattering 

= 3 if total, scattering, and y*fission. 

NINC 

Index for locating point cross sections. Used as flag in 
FISGEN to determine if transport process has begun. 
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TABLE XI. 


Definitions of Variables in Common PDET 

Variable Definition 

ND Number of detectors (2 1). 

NNE Number of neutron energy bins. 

NE Total number of energy bins. 

NT Number of time bins. 

NA Number of angle bins. 

NRESP Number of energy response functions 1). 

NEX Number of extra arrays of length NMTG to be set aside. 

NEXND Number of extra arrays of length ND to be set aside. 

NEND NE*ND. 

NDNR ND*NRESP. 

NTNR NT*NRESP. 

NTNE NT*NE. 

NANE NA*NE. 

NTNDNR NT*ND*NRESP. 

NTNEND NT*NE*ND. 

NANEND NA*NE*ND. 

L0CRSP Location of cell zero of response functions. 

L0CXD Location of cell zero of detector positions. 

L0CIB Location of cell zero of energy bins. 

L0CC0 Location of cell zero of angle (cosine) bins. 

L0CT Location of cell zero of time bins. 

L0CUD Location of cell zero of array UD. 

L0CSD Location of cell zero of array SD. 

L0CQE Location of cell zero of array QE. 

L0CQT Location of cell zero of array QT. 
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TABLE XI. (Cont'd.) 


Variable 


Definition 


L0CQTE Location of cell zero of array QTE. 

L0CQAE Location of cell zero of array QAE. 

LMAX Last cell used in blank common. 

EFIRST Upper energy limit of first energy bin. 

EGT0P Upper energy limit of first gamma-ray bin. 
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TABLE Xn. 


Definitions of Variables in Common FISBNK 


Variable Definition 


MFISTP Index to indicate that fissions are to be treated if > 0. 

NFISBM Location in blank common of cell zero of the fission 

bank. 

NFISH Number of fissions produced in the previous batch 

(generation). 

FT0TL Total fission weight from all collisions in this batch 

(generation). 

FWATE Total weight of fission neutrons stored in bank in this 

batch (generation). 

WATEF Weight of fission neutron stored in bank. 







TABLE Xm. 


L-efinitions of Variables in Common GE0MC 
as foimd in 05R Geometry 


Variable 


Definition 


X2, Y2, Z2 

XI, Yl, Z1 
ETA 

ETAUSD 

IBLZ 

IBZN 

MARK 


NMED 

NREG 


Coordinates at tentative end-of-flight or if the trajectory 
is in an internal void; X2, Y2, Z2 are the direction 
cosines of the trajectory. 

Starting coordinates for the particle. 

Number of mean free paths to be traversed if flight goes 
to X2, Y2, Z2. 

Number of mean free paths actually traver sed after the 
call to GE0M. 

An index to the medium number for the special geometry 
packages. For GENERAL GE0M, IBLZ is a packed 
work giving the block and zone of the end of flight. 

A dummy variable. 

A flag set by GE0M indicating the results of the trajec¬ 
tory calculation 

= 1 for completed flight 
= 0 for boundary crossing 
= -1 for escape 

= -2 for entering an internal void. 

Medium number at end of flight or of medium about to 
be entered at a boundary crossing. 

Region number at end of flight; not set at boundary 
crossings. 






TABLE XIV. 


Definitions of Variables in Common PAREM 
as found in Combinatorial Geometry 

Variable Definition 

XB(3) Coordinates of the starting point of the present path. 

WB(3) Direction cosines of particle trajectory. Equal to U, 

V, andW. 

IR Combinatorial zone of present particle position. 

WP(3) Temporary storage of WB(3). 

XP(31 Temporary storage of XB(3). 

IDBG Set non-zero to initialize a debug printout. 

IRPRIM Next region to be entered after a call of Gl. 

NASC Body number of last calculated intersection. Set negative 

to indicate source or collision point not on a body sur¬ 
face. 

LSURF Surface of body NASC where intersection occurred. 

Positive if particle is entering the body and negative 
when exiting. 

NB0 Body number and a sign used to define zones. Input in 

zone description as positive when zone is contained in 
body and as negative if zone is outside body. 

LRI Entry surface calculated in GG. 

LR0 Exit surface calculated in GG. 

RIN Distance to entry calculated in GG. 

R0UT Distance to exit calculated in GG. 

KL00P Trajectory index of present path incremented in Gl. 

L00P Index of last trajectory calculated for body NB0. If 

L00P is equal to KL00P, GG returns immediately with 
old values of RIN, R0UT, LRI, and LR0. 

Body type of body NB0 (indicates B0X, SPH, etc.). 


ITYPE 
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TABLE XIV. (Cont’d.) 

Variable Definition 


PINF Machine infinity. 

N0A Not used. 

DISV Distance from XB(3) to present point. 












TABLE XV. 

* 

Definitions of Variables in Common 0RGI 

Variable Definition 

M.‘L. jC Set 1 in G1 if trajectory end point is reached before 

next intersection. Otherwise set to 0. 

DISTO Distance from point XB(3) to next scattering point. 

Used in G1 to avoid calculating the next zone if a 
scattering event occurs before the intersection. 

NMEDG Zone number IR from a L00KZ call. Stored in BLZNT 

by MS0UR. 

_ 

Note: Variabl names are not the same in all routines. 
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TABLE XVI. 


Definitions of Random Walk Variables in Blank Common 


Mnemonic 

Variable Definition 

Name 


ENER(IG) 

VEL(IG) 


FS(IG) 


BFS(IG) 

WTHI(IG, 

NREG) 

WTL0(IG, 

NREG) 

WTAV(IG, 

NREG) 

PATH(IG, 

NREG) 

NSPL(IG, 

NREG) 

WSPL(IG, 

NREG) 

N0SP(IG, 

NREG) 

WN0S(IG, 

NREG) 

RRKL(IG, 

NREG) 


Upper energy boundary of group IG (in eV). 

Velocity corresponding to the mean energy for neutron 
groups and the speed of light for gamma-ray groups 
(in cm/sec). 

Unbiased source spectrum—unnormalized fraction of 
source particles in each energy group--transformed to 
c.d. f. by S0RIN. 

Biased source spectrum—relative importance of each 
energy group—transformed to biased c. d. f. by S0RIN. 

Weight above which splitting is performed (vs. group 
and region). 

Weight below which Russian roulette is performed (vs. 
group and region). 

Weight tc be assigned Russian roulette survivors (vs. 
group and region). 

Exponential transform parameters (vs. group and region). 
Splitting counter (vs. group and region). 

Weight equivalent to NSPL. 

Counter for full bank when splitting was requested (vs. 
group and region). 

Weight equivalent to N0SP. 

Russian roulette death counter (vs. group and region). 
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TABLE XVI. (Cont'd.) 


Mnemonic 

Variable Definition 

Name 


WRKL(IG, 

NREG) Weight equivalent to RRKL. 

RRSU(IG, Russian roulette survival counter (vs. group and 

NREG) region). 

WRSU(IG, 

NREG) Weight equivalent to RRSU. 

INIWHI (IG, 

NREG) Initial values of WTHI array. 

INIWL0(IG, 

NREG) Initial values of WTL0 array. 

INIWAV(IG, 

NREG) Initial values of WTAV array. 

FWL0(NREG) Weights to be assigned to fission daughters (vs. 

region). 

INIFL0(NREG) Initial values of FWL0. 

GWL0(IG, Weights to be assigned to secondary particles (vs. 

NREG) group and region). 

EPRB(IG, Relative importance of energy groups after scatter- 

NREG) ing (vs. group and region). 

NSCT(IG, 

NREG) Number of real scatterings (vs. group and region). 

WSCT(IG, 

NREG) Weight equivalent to NSCT. 

NALB(IG, Number of albedo scatterings (vs. group and 

NREG) region). 

WALB(IG, 

NREG) Weight equivalent to NALB. 

NFIZ(IG, 

NREG) Number of fissions (vs. group and region). 





tABLfifcVI. (Cont'd.) 


Mnemonic » 

Variable Definition 

Name 


WFIZ(IG, 

NREG) 

NGAM(IG, 

NREG) 

WGAM(IG S 
,NREG) 4 

NSCA(IMED) 

FISH(IG, 

IMED) 

FSE(IG, 

IMkD) 

GMGN(IG, 

IMED) 


Weight equivalent to NFIZ. 

Number of secondary productions (vs. group and region). 

Weight equivalent to NGAM. ' 

Scattering counter (vs. cross-section medium). 

Probability of generating fission neutron (vs. group 
and medium). 

Source spectrum for fis$ion-induced neutrons for each 
group—input as frequency of group IG. 

Probability of generating secondary particle (vs. group 
and medium). 


\ 
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IV. THE PICTURE PROGRAM 

\ 

This program is used to obtain printer plots of two-dimensional 
slices through the three-dimensional geometry described by the combi¬ 
natorial geometry (CG) package. The 05R geometry version^°| has 
been modified for use with CG and some minor changes to the program 
have been made. It is now possible to obtain a picture the full width of 
a computer page and the logic for determining the medium and/or region 
number has been changed so that a test of the geometry ray tracing cap¬ 
ability is also made. That is. instead of just testing the routine L00KZ, 
this modification forces the geometry package to execute the routines 
and logic used during the random walk process. This program has been 
very useful in debugging geometry input. It is suggested tjiat any com¬ 
plex geometry be tried in this program before using it in the M0RSE 
code. A sample problem is provided. Besides the geometry package, 
a main program, PICTURE, and two routines, MESH and PRINT, are re¬ 
quired. A discussion of these routines is included and a definition of 
the variables in common, PICT, is given in Table XVII. 


Preceding page blank 
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4.1 Program PICTURE 

PICTURE is the executive routine for the PICTURE program 
and reads in the input, calculates the coordinates of the picture to be 
plotted and controls the calls to other routines. There are several dif¬ 
ferent ways in which a two-dimensional slice through the geometry may 
be obtained. These different options are discussed in Appendix B with 
discussions of the input. The characters to be printed for corresponding 
media or regions may be changed by altering the values in ATABLE as 
given in the data statement. 

Subroutines called: J0MIN, PRINT 
Commons required: PICT 

Variables required: Several input cards are read. 

Variables changed: All variables in common PICT. 

Significant internal variables: 

INT - input logical unit, 

I0T - output logical unit, 

NADD(l) - first location in blank common for storage 
of geometry data, . 

flag to indicate which set of input options 
was used to define the two-dimensional 
slice. 


NCK 






Program PICTURE 


START 


[Set input/output devices to 5 an d 6. 
j 1/0: Sym bol optionT] 

| Read symbol table if required 1 


Call J0MIN to input geometry data.] 


1/0: Options, number of pictures for this 
geometry, media or region title. _ 

1/0: picture coordinates upper left and 
_ lower right corner. _ 

fl/0: direction numbers for U and v nve<n 


1 ^/ 0 - interval number and 


Were both upper left and lower\„ - 

right coordinates (riven? 1 * es —* Calculate 

-- J DELU given 


No 

,KCK=1 


Use i nput DELU or DELV values and NU. NvT) 


NU cr NV 
given DELU 
I repeat for V 


c. !! PRINT 


f Is there ano ther geometry'input?] 


3 










4. 2 Su broutine PRINT (KXX, KYY, ATABLE) 

This routine controls the printing of the picture. First it is 
determined if the picture is to be more than one "page” (130 columns) 
wide and then for each line of the picture on the first page the infor¬ 
mation to be printed is determined by calling MESH and printed with a 
130A1 format. The next lines are then calculated and printed until that 
"page” is finished. Note that a page refers to width not length. Thus, 
as much detail may be obtained as necessary in both directions by piec¬ 
ing together the output. 

Called from: PICTURE 

Subroutine called: MESH 
Commons required: PICT 
Variables required: 

KXX - number of intervals in the U direction 

(direction of paper movement through 
the printer), 

KYY - number of intervals in the V direction 

(line), 

ATABLE - table of characters to be printed. 

Significant interval variable: 

NPAGES - number of subpictures required to cover 
the width of the total picture, 

I0T - output logical unit, 

NV - number of characters per line (characters/ 

page width). 
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*1 


Subroutine PRINT 



RETURN 









4.3 


Subroutine MESH (XS, YS, ZS, NV) 


Subroutine MESH is used by the PICTURE package to set up 
one line of print in the array NST0R. Both for efficiency and to debug 
the combinatorial geometry package, this version has been modified to 
work exactly like particle tracking. L00KZ is first called to deter¬ 
mine the zone of the first grid point. A trajectory to the last grid point 
is then initialized, and successive calls to G1 ’’track” a particle to the 
last point, setting the region of each grid point in NST0R. By setting 
IRG negative, o ro, or positive, either NREG, IR, or NMED will be 
stored in the print array NST0R. 


Called from: PRINT 

Subroutines called: L00KZ, G1 

Commons required: PICT, G0ML0C, PAREM, 0RGI, blank. 
Variables required: 


XS, YS, ZS 
NV 
DELV 
IRG 


coordinates of first grid point, 

number of grid points, 

distance between grid points, 

flag to print NREG, IR, or NMED if IRG 
is negative, zero, or positive. . 


Variables changed: NST0R - print array. 


Significant internal variables: 

J - grid point index, 

IST0R - value to be stored in NST0R between 

successive boundary crossings. 
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Subroutine MESH 



No 



_ 

Loop from next (c 

— 
last .grid points. 


| Store IST0R in NST0R array. | 


_ 


|Set next point t 

3 last point +1. I 


fj R » IRPRIm] 


(Is traj ect ory complete? ) -No- 

i Yes 


i 

! 


RETUR! 





TABLE XVH. 


* 

Definition of Variables in Common PICT 


Variable 

Definition 

DELU 

The increment in geometry units between lines in the 
picture in the U direction. 

DELV 

The increment in geometry units between lines in the 
picture in the V direction. 

XO, YO, ZO 

The coordinates in geometry units defining the first 
point (upper left hand corner) of the picture. 

XU, YU, ZU 

The length (in direction U) of the picture in geometry 
units. 

XV, YV, ZV 

The width (in direction V) of the picture in geometry 
units. 

NST0R(13O) 

An array used to store the medium, region, or zone 
number for one line of the picture. 

IRG 

A flag indicating that region, zone, or medium geo¬ 
metry parameter should be printed if IRG is negative, 
zero, or positive, respectively. 


Used in the PICTURE program. 
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V. M0RSE SAMPLE PROBLEMS 

5.1 Sample Problem Number 1 

(7) 

This sample problem is similar to that used previously 
which calculates neutrons and secondary gamma rays in infinite homo¬ 
geneous air. The source is a 12. 2- to 15-MeV neutron source. Two 
types of estimators are used in the analysis of the particle histories. 
First, a boundary-crossing estimator is used for five spherical shells. 
(This part of the problem is similar to the M0RSE sample problem.) 
Second, a collision density estimator is used in RELC0L for nine de¬ 
tector regions (detectors 6-14). By choosing the spherical shells as 
the midpoints of the regions, an internal check is provided. Also, by 
choosing the geometry regions to correspond to the collision density 
detectors, the results from SAMB0 can be compared with results from 
ENDRUN. With the inclusion of the combinatorial geometry package, 
an additional estimator based on track length/unit volume has been 
added with results also output in ENDRUN. Both the collision density 
and track length per unit volume estimators are obtained for the geo¬ 
metry regions without s. ecific coding. 

Since the collision density estimator in RELC0L scores 

2 

4ffR. <p there will be some differences between these estimates and 

those from the collision ..ensity estimator output in ENDRUN which 
2 

scores 4ir <R> <p. (The symbols < > denote an average over a spatial 

region.) 

Table XVIII illustrates the comparison for one of several 
sets of answers as calculated by the various estimators. Descriptions 
of the user routines INSC0R, RELC0L, and BDRYX required for the 
sample problem are given. 
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TABLE XVm. 


Comparison of Results of Sample Problem 1 


Detector No. 

Type and Location 
of Detector 

Result 

Response #1 Group 1 * 

4ttr 2 fluence fluence/eV 
top group 

3 

Boundary crossing 
at 450 meters. 

8.19(-2) 

2. 92(-8)*A 
= 8.19(-2) 

9 

Collision density for 
volume between 400- 
500 meters. 

7. 59(-2) 

2. 71(-8)*A 
= 7. 59(-2) 

Region 4 

Volume between 
400-500 meters 
(ENDRUN results). 

Real 

Collision 
Region 4 

Top Group 

1.11(-18)*A*B 
= 7. 76(-2) 


Volume between 
400-500 meters 
(ENDRUN results). 

Track 

Length 

Region 4 

Top Group 

1.18(-18)*A*B 
= 8. 36(-2) 


A = width of the top group 2. 8(+6) eV. 


B = 4ff (R: ^ for 4^ region = 2. 5(+10). 






5.1.1 Subroutine INSC0R - Sample Problem Number 1 

This user routine is used to calculate the volume of each of 
the detector regions (detectors 6-14) in the sample problem. The re¬ 
ciprocal of the volume is stored in the array in blank common set aside 
for FACT, a detector-dependent normalization. Because two types of 
estimators are used in the problem, normalization values of 1 are 
stored for the boundary detectors. The detector number NDC for the 
last boundary crossing detector (in this case, 5) must be input. 

Called from: SC0RIN 

Subroutines called: None 

Commons required: PDET, USER, DETCUT, blank. 

Variables required: 

NDC - the number of the last boundary crossing 

detector. Detector radii from blank 
common. 

Variables changed: 

Normalization array in blank common. 


Ill 










Subroutine INSC0R 
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5.1.2 Subroutine RELC0L - Sample Problem Number 1 

This subroutine estimates the iluence from the density of 

2 

collisions. For comparison purposes, 4trR times the fluence is stored 
where R. is the radius to each collision site. Thus, the contribution is 

4t p. 2 

1 £ m 


where 


R. 


-- X 2 + Y 2 + 

3 


z 2 , 

1 1 ’ 


WTBC = particle weight before collision, 

Is = cross section for the group (IG0) of the incoming 
particle. 


A detector number from NDC + 1 to ND is chosen depending 
on the value of R i . 

Called from: BANKR(5) 

Subroutines called: FLUXST, NSIGTA 

Commons required: PDET, NUTR0N, USER, DETCUT 

Variables required: 

X, Y, Z, IG0, NMED, AGE, WTBC from NUTR0N common 
TS1G and detector radii. 

Significant internal variables:- 

R - radial distance to collision site, 

C0N - fluence estimate. 


113 








l 

Subroutine RELC0L 
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5.1. 3 Subroutine BDRYX - Sample Problem Number 1 ! 

I 

This subroutine is called whenever the partible in the random 
walk encounters a change in geometry media. If the source-to-collision 
distance corresponds to a detector position, the reciprocal of the cosine 
of the angle f^om the radius vector is used as a fluence estimate. 


Called from: 

BANKR(7) 

Subroutines called: 

ERR0R 

Coriimons required: 

Variables required: 

USER, 

NUTR0N, DET 

X, Y, Z, 

u, V, 

W, WATE (from common NUTR0N). 

Significant internal variables: 

R21 

- 

iradial distance to boundary crossing, 

' i 

R2 

- 

99% of R21, 

R22 

- 

101% of R21. 

C0S 

- 

cosine of angle between particle direction 
and radius vector, 

ABC0S 

- 

absolute value of C0S, 

C0N 

- 

fluence estimate. 


i 


\ 

1 


i 








Subroutine BDRYX 
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5.2 


Sample Problem Number 2 


This sample problem is used to illustrate the next event 

estimator. The neutron and secondary gamma ray spectra at several 

4 

point detectors located 10 cm from a polyethylene slab is calculated. 
The source is a beam of neutrons with a fission spectrum energy dis¬ 
tribution. The use of the point detector routines to estimate the energy 
dependent fluence is illustrated and in particular, for secondary gamma 
rays the contribution of the uncollided or first flight gamma rays are 
scored in subroutine SGAM. The estimates from real collisions are 

made in subroutine RELC0L. Results are recorded in Table XIX in 
2 

units of particles cm source neutron. 


TABLE XIX. 

Results of Last Flight Estimator 
Sample Problem 


Detector 

2 

Neutrons/cm 

Per 

Source Neutron 

2 

Gammas/cm 

Per 

Source Neutron 

1 

6. 24(-l1) 

2. 36(-ll) 

2 

6. 44(-10) 

4. 82(-ll) 

3 

1. 32(-9) 

4. 93(-ll) 

4 

1. 56(-9) 

4. 94(-ll) 

5 

6. 92(-l) 

2.15(-3) 
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5. 3 Sample Problem Number 3 

This sample problem serves two purposes by (1) demonstrat¬ 
ing the use of the PICTURE code, and (2) verifying the tracking capability ■ 
of the combinatorial geometry package in a moderately complex three- 
dimensional system. In the PICTURE program, subroutine G1 is called 
from MESH in the exact manner that it is called from G0MST during a 
M0RSE calculation. As such, PICTURE works by tracking rays and 
thus can be used to debug a combinatorial geometry input. 

The sample problem itself is a picture of a tank. The tank 
model itself was constructed purely as an illustration of the combina¬ 
torial geometry and is in no way accurate or detailed. Notice that the 
use of both the 0R operator and the ARB body is demonstrated. 
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Output from PICTURE Sample Problem 
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TABLE XX. (Cont’d.) 
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APPENDIX A 


ENERGY INDEXING SCHEME 

For most problems the user does not need to worry about the 
energy-indexing scheme used in the M0RSE Monte Carlo code. Energies 
corresponding to the group boundaries in the forward group structure 
are input and the code takes care of the rest. The analysis package 
SAMB0 performs the proper bookkeeping for labeling. A diagram for 
the energy group structure is given in Fig. Al. However, some clarifi¬ 
cation of the indexing scheme for several types of problems is advanta¬ 
geous. This clarification can perhaps best be made through the use of 
examples. Consider a five-group coupled neutron-gamma-ray problem 
with the following group boundaries: 


Group 

Energy (eV) 

with NMGP = 3, NMTG = 5, 

1 

15(+6) 

NGPQTN = 3, NGPQTG = 2 

2 

5(+6) 


3 

l(+6) 


EB0TN 

3 (+3) 


4 

10 (+6) 


5 

2 (+6) 


EB0TG 

4(+5) 



There are several options which might be considered; namely, 
(1) a forward coupled problem, (2) a forward neutron-only problem, (3) 
a forward gamma-ray-only problem, (4) an adjoint coupled problem, (5) 
an adjoint neutron-only problem, and (6) an adjoint gamma-ray-only 
problem. Table AI gives the energies as they are indexed during a 
M0RSE run. Table All gives the values of variables NQTl, NQT2, and 
NQT3 as they appear in USER common for the same six problems. 

Preceding page blank 
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Forward 


Adjoint 


NGPQT1 


NGPQT2 


NGPQT3 


1 


NMTG 



- 

NGPQTN 



- 


NMTG - NGPQTN 

- 


- 

NMGP 






- 


NMTG - NMGP 

- 


- 

NMGP + NGPQTG 



- 


NMTG - NMGP - NGPQTG 



- 

NMTG 


1 


Fig. Al. Diagram of energy group structure. 


NGPQT3 


NGPQT2 


NGPQT1 


-NGPQTO 


NGPQT3 = NGPQT1 for neutron only or gamma ray only problem. 
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APPENDIX B 
INPUT INSTRUCTIONS 


B. 1 M0RSE Input Instructions 

The input read 1 by subroutine INPUT is as follows: 

CARD A (20A4) 

Title card. (Any character other than a blank or alphameric in 
column one will terminate the job.) 


CARD B (1315) 
NSTRT 
NM0ST 


NITS 

NQUIT 

NGPQTN* 

NGPQTG* 

NMGP* 


number of particles per batch, 

maximum number of particles allowed for in the \ 
bank(s); may equal NSTRT if there is no splitting, 
fission, and secondary generation during ex¬ 
ecution. If bank size is exceeded by more than 
50 due to fission or secondary gamma ray gen¬ 
eration the job is terminated, 

number of batches, 

number of sets of NITS batches to be run without 
calling subroutine INPUT, 

number of neutron groups being analyzed, 

number of gamma-ray groups being analyzed, 

number of primary particle groups for which 
cross sections are stored; should be same as 
NGP (or the same as NGG when NGP = 0) on 
Card XB read by subroutine XSEC, 


Input and output logical unit members are defined in MAIN and ii Gl. 
Unit 16 is a temporary storage unit used by J0M3N. 

See Table BI for sample input. 
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* I 

NMTG - total number of groups for which cross sections 
are stored; should be same $s NGP + NGG as 
read on Card XB read by subroutine XSEC, 

NC0LTP 4 set greater than zero if a collision tape is de¬ 
sired; the collision tape- is written by the user 
routine BANKR, 

IADJM - set greater than zero for an adjoint problem, 

MAXTiM - maximum clock time in minutes allowed for the 

I problem to be on the computer, 

]MEDIA - number of cross-section media; should agree | 
v'th NMED cn Card XB read by subroutine XSEC, 

MEDALB - albedo scattering medium is absolute value oi 
MEDALB; if 

1 = 0, no albedo information to be read in, 

<0, albedo only problem—no cross sections 
are to be read, 

> 0, coupled albedo and transport problem. 


CARD C (415, 5E10. 5) 


IS0UR - source energy group if > 0, if IS0UR £ 0, 

S0RIN is called for input of Cards El and E2, 

NGPFS - rumlber of groups for which the source spectrum 
is to be defined. If IS0UR £ 0, NGPFS * 2, 

ISBIAS - \ no source energy biasing if set equal to zero; 

otherwise the source energy is to be biased, 
and Cards E2 are required, 


N0TUSD 

WTSTRT 


an unused variable, 

! 

weight assigned to each source particle, 


EB0TN 


lower energy limit of lowest neutron group (eV) 
(group NMGP), 


EB0TG 


lower energy limit of lowest gamma-ray group 
(eV) (group NMTG), : 


\ 


* 

See Table BI for sample input; 
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TCUT - age in sec at which particles are retired; if 
TCUT = 0, no time kill is performed, 

VELTH - velocity of group NMGP when NGPQTN > 0; i. e., 
thermal-neutron velocity (cm/sec). 


CARD D (7E10.4) 

XSTRT ) 1 

YSTRT > - coordinates for source particles, 
ZSTRT ' 


^GSTRT 

UINP ) 
VINP > 
WINP ) 


starting age for source particles, 

source particle direction cosines; if all are 
zero, isotropic directions are chosen. 


\ 

Source data on Cards C and D, may be overridden by changes 
in subroutine S0 URCeL 


CARDS El (7E10.4) (omit if IS0UR on Card C > 0) 

NGPFS values of FS, where FS equal the unnormalized fraction 
of source particles in each group. 


CARDS E2 (7E10. 4) (omit if IS0UR > 0 or if IS0UR * 0 and LSBIAS = 0) 

If ISBLAS >0, NGPFS values of BFS, the relative importance of 
a source in group I, are required. 

CARDS F (7E10. 4) 

NMTG values of ENER, the energies (in eV) at the upper edge of 
the energy group boundaries. (Note: The lower energies of 
groups NMGP and NMTG were read on Card C.) 


CARD G (215, 5X, 3611, 1311) (omit if NC0LTP on Card B * 0) 

NHISTR - logical tape number for the first collision tape, 

NHISMX - the highest logical number that a collision tape 

I may be assigned, 
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NBIND(J), J-l, 36 - an index to indicate the collision para¬ 

meters to be written on tape, 

NC0LLS(J), J=l, 13 - an index to indicate the types of collisions 
to be put on tape. (See Table BII for in¬ 
formation concerning NBIND’s and 
NC0LLS.) 


CARD H (Z12) 

RAND0M - starting random number. 


CARD I (715) 
NSPLT 
NKILL 

NPAST 

N0LEAK 

IEBIAS 

MXREG 1 

MAXGP 


index indicating that splitting is allowed if > 0, 

index indicating that Russian roulette is allowed 
if >0, 

index indicating that exponential transform is in 
voked if > 0 (subroutine DIREC required), 

index indicating that non-leakage is invoked if 

> 0 , 

index indicating that energy biasing is allowed 
if >0, 

number of regions described by geometry input 
(’•ill be set to one if * 0), 

group number of last group for which Russian 
roulette or splitting or exponential transform is 
to be performed. 


I' ENDRUN is used, a data array relating media numbers to region 
numbers must be given in a data statement in ENDRUN. 
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'J, 


1 


CARD J (615, 4E10. 5) (omit if NSPLT + NKILL + NPAST = 0) 

from energy group NGPl to energy group NGP2, 
inclusive, in steps of NDG and from region NRG1 
to NRG2, inclusive, in steps of NDRG, the follow¬ 
ing weight standards and path-stretching para¬ 
meters are assigned. If NGPl = 0, groups 1 to 
MAXGP will be used; if NRG1 = 0, regions 1 to 
MXREG will be used (both in steps of one). Us¬ 
ually NDG = 1 and NDRG = 1. 

weight above which splitting will occur, 

weight below which Russian roulette is played, 

weight given those particles surviving Russian 
roulette, 

path-length stretching par. meters for use in ex ¬ 
ponential transform (usually 0 £ PATH < 1). 

The above information is repeated until data for all groups and 
regions are input. 

End Cards J with negative value of NGPl (ex., -1 in columns 4 
and 5). 

CARDS K (7E10.4) (omit if IEBIAS on Card I s 0) 

((EPR0B(IG, NREG), IG = 1, NMTG), NREG = 1, MXREG) 

Values of the relative energy importance of particles leaving a 
collision in region NREG. Input for each region must start on a 
new card. 

CARD L (415) 

NS0UR - set £ 0 for a fixed source problem; otherwise the 
source is from fissions generated in a previous 
batch, 

MFISTP - index for fission problem, if £ 0 no fissions are 
allowed, 

NKCALC - the number of the first batch to be included in the 
estimate of k; if * 0 no estimate of k is made, 


NGPl 

NDG 

NGP2 

NRG1 

NDRG 

NRG2 


WTHIH 

WTL0W 

WTAVE 


PATH 
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N0RMF 


the weight standards and fission weights are un¬ 
changed if £ 0; otherwise fission weights will be 
multiplied, at the end of each batch, by the latest 
estimate of k and the weight standards are mul¬ 
tiplied by the ratio of fission weights produced in 
previous batch to the average starting weight for 
the previous batch. For time-dependent decay¬ 
ing systems, N0RMF should be > 0. 


CARDS M (7E10.4) (omit if MFISTP on Card L * 0) 

(FWL0(I), 1=1, MXREG) values of the weight to be assigned to 
fission neutrons. 


CARDS N (7E10.4) (omit if MFISTP on Card L^O) 

(FSE(IG, IMED), IG = 1, NMGP), IMED = 1, MEDIA) the fraction 
of fission-induced source particles in group IG and medium IMED. 

Note: Input for each medium must start on a new card. 


CARDS 0 (7E10. 5) (omit if NGPQTN = 0 or NGPQTG = 0, i. e., include 
if coupled neutron-gamma-ray problem) 

((GWL0(IG, NREG), IG = 1, NMGP or NMTG - NMGP), NREG = 1, 

MXRE< ^ _ average number of secondary particles to be 
generated. NMGP groups are read for each 
region in a forward problem and NMYG - 
NMGP for an adjoint. Input for each region 
must start on a new card. 
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TABLE BI. 


Sample Group Input Numbers for Some 
Representative Problems* 


Case A - Neutron Only Cross Sections (22 groups) 

Case B - Gamma-Ray Onlv Cross Sections (18 groups) 

Case C - Neutron-Gamma-Ray-Coupled Cross Sections (22-18 groups) 


Input Variable 

Problem Type 

Case A 

Top 14 groups 

Case B 

Top 17 groups 

Case C 

Neutrons Only 

Top 14 groups 

Case C 

Gamma Rays Only 
Top 17 groups 

Case C 

Neutron-Gamma 
Top 14 Neutron, 
Top 17 Gamma 
groups 


M0RSE Input: 







NGPQTN 

14 

0 

14 

0 

14 


NGPQTG 

0 

17 

0 

17 

17 

CARD B, 

NMGP 

22 

18 

22 

18 

22 

Variables 

NMTG 

22 

18 

22 

18 

40 


NGP 

22 

18 

22 

0 

22 


NGG 

0 

0 

0 

18 

18 

CARD XB 

INGP 

22 

18 

40 

40 

40 

Variables 


For cross sections with full downscatter, NDS = NGP, NDSG = NGG, 
INDS = iNGP, and ITBL = number of downscatters + number of up- 
scatters + 3. Usually, ISGG = number of upscatters + 4. 
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TABLE BE. 

BANKR Arguments (NC0LLS) 


BANKR 

Argument 

Called 

from 

Location of Call in Walk 

-1 

M0RSE 

After call to INPUT—to set parameters 
for new problem 

-2 

M0RSE 

At the beginning of each batch of NSTRT 
particles 

-3 

M0RSE 

At the end of each batch of NSTRT 
particles 

-4 

M0RSE 

At the end of each set of NITS batches— 
a new problem is about to begin 

1 

MS0UR 

After a source event 

2 

TESTW 

After a splitting has occurred—com¬ 
mented in column 1 

3 

FPR0B 

After a fission has occurred 

4 

GST0RE 

After a secondary particle has been 
generated 

5 

M0RSE 

After a real collision has occurred— 
post-collision parameters are avail¬ 
able 

6 

M0RSE 

After an albedo collision has occurred— 
post-collision parameters are available 

7 

NXTC0L 

After a boundary crossing occurs (the 
track has encountered a new geometry 
medium other than the albedo or void 
media) 

8 

NXTC0L 

After an escape occurs (the geometry 
has encountered medium zero) 

9 

M0RSE 

After the post-collision energy group 
exceeds the maximum desired—com¬ 
mented in column 1 
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TABLE BH. (Cont'd.) 


BANKR 

Argument 

Called 

from 

Location of Call in Walk 

10 

M0RSE 

After the maximum chronological age 
has been exceeded—commented in 
column 1 

11 

TESTW 

After a Russian roulette kill occurs— 
commented in column 1 

12 

TESTW 

After a Russian roulette survival 
occurs—commented in column 1 

13 

GST0RE 

After a secondary particle has been 
generated but no room in the bank is 
available—commented in column 1. 
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B. 2 05R Geometry Input Instructions 

Read by subroutine J0MIN and the specific input depends on the 
geometry packages used. 


B. 2.1 Spherical GE0M 
CARD GA (15, DIO. 5) 

MED - medium number interior to R(> 0), 

R - outer radius of sphere or spherical shell contain¬ 

ing MED. 

Repeat Card GA for all radii (s 20) in increasing order. End 
Card GA input with blank card. 


CARD GB (DIO. 5) 

R - region radius of sphere or spherical shell con¬ 

taining regions. Region numbers are assigned 
in consecutive order starting with 1, and R must 
be in increasing order. 

Repeat Card GB for all radii (s 20) in increasing order. End Card 
GB input with blank card. If no regions are desired, a blank card 
must be used to signal no region geometry. 

B.2.2 Slab GE0M 

CARD GA (15, DIO. 5) 

MED - medium with Z as lower bound (> 0), 

Z - lower limit of medium MED. 

Repeat Card GA for all boundaries with the last card containing 
MED = 0 and the boundary of the system. 







CARD GB (DIO. 5) 

Z - lower limit of region boundary. Region numbers 

are assigned in consecutive order starting with 
1, and Z must be in increasing order. 

Repeat Card GB for all region boundaries. 

End Card GB input with a blank card. If no region geometry is 
desired, a blank card is required. 


CARD GC (4D10. 5) 


XL 

lower 

XU 

upper 

YL 

lower 

YU 

upper 


boundary of system in X direction, 
boundary of system in X direction, 
boundary of system in Y direction, 
boundary of system in Y direction. 


B. 2.3 Cylindrical GE0M 
CARD GA (15, 5X, A8) 

NREGIN - flag to indicate material media (=2) or both 
region and material media (=1), 

SEX - sex of programmer. 

CARD GB (E10. 5) 

R - radii of the cylindrical shells describing the 

material media in ascending order. 

Repeat Card GB until all radii have been input. 

End Card GB input with a blank card. 
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CARD GC (ElO. 5, 1215/815) 


H - upper height of medium M(I) (>0), cylinders as¬ 

sumed to start at H = 0, 

M(I) - media for the cylindrical shells for this height. 

Repeat Card GC until all height intervals have been input. 

End Card GC input with a blank card or if there are more than 12 
radial intervals, 2 blank cards. 


CARD GC (ElO. 5) (omit if NREGIN = 2) 

RG - radii of the cylindrical shells describing the 

region geometry in ascending order. 

Repeat Card GD until all region geometry has been input. 

End Card GD with a blank card. 


CARD GE (ElO. 5, 1215/815) (omit if NREGIN = 2) 

HG - upper height of region MG(I), 

MG(I) - region numbers for the cylindrical shells for this 

height. 

Repeat Card GE until all height intervals have been input. 

End Card GE input with a blank card or if there are more than 12 
radial intervals, 2 blank cards. 


B. 2.4 General GE0M 


CARD GA (15, 5X, A6, 
NSTAT 

SEX 

STATUS 


IX, A7) hollerith left adjusted 

flag to indicate material media only if 2 and both 
region and material media if 1, 

sex of the programmer (select one from MALE, 
FEMALE, or blank indicating uncertain), 

marital status of programmer. 
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CARD GB (2A4, A3, 5(D10. 5, Al)) 

DUMMY(3) - bollerith characters not used, 

FIN(I) - zone boundaries increasing order along the x 
axis, 

BCD(I) - flag to indicate end of input if blank, comma 
means to continue. 

Repeat in format (6(D10. 5, Al))if more than five boundaries 
along the x axis are needed. 

CARD GC - same as CARD GB except for y axis. 

Repeat in format (6(DIO. 5, Al)) if more than five boundaries 
along the y axis are needed. 

CARD GD - same as Card GB except for z axis. 

Repeat in format (6(D10. 5, Al)) if more than five boundaries 
along the z axis are needed. 

CARD GE (A4, A2, 315) 

BCD1 - hollerith Z0NE, 

BCD2 - dummy, 

integers which specify the zone as being the 
NXZN0th zone in the x direction, NYZN0th zone 
in the y direction, and NZZN0th zone in the z 
direction. 

CARD GF (2A4, A3, 5(D10. 5, Al)) 

DUMMY(3) - hollerith characters not used. 

FIN(I) - block boundaries in increasing order along the 


NXZN0 

NYXN0 

NZZN0 









BCD(I) - flag to indicate end of input if blank, comma 
means to continue. 

Repeat in format (6(D10. 5, Al)) if more than five boundaries 
along the x axis are needed. 


CARD GG - same as Card GF except for y axis. 

Repeat in format (6(D10. 5, Al)) if more than five boundaries 
along the y axis are needed. 

CARD GH - same as Card GF except for z axis. 

Repeat in 7 ormat (6(DIO. 5, Al)) if more than five boundaries 
along the z axis are needed. 


CARDS GI to G0 describe the geometry for a block and must be included 
for each block in me zone. 


CARD GI (A4, A2, 315) 


BCD1 

- 

holleriih BL0C, 

BCD2 

- 

dummy, 

NXBND 

NYBND 

NZBND 

1 - 

integers which specify the block as being the 
NXBNDth in the x direction, the NYBND^h in the 
y direction, and the NZBND in the z direction. 

CARD GJ (3A4, 

10(15, 

Al)) 

NAM2 

- 

hollerith MEDI, 

DUM(2) 

- 

dummy, 

INP(I) 

- 

a list of media sector by sector in the block, 

BCD(I) 

- 

flag to indicate end of input if blank, a comma 
means to continue. 


Continuation with 12(15, Al). 
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CARD GK (3A4, 10(15, Al)) 


NAM2 - hollerith SURF, 

DUM(2) - dummy, 

INP(I) - a list of quadratic surfaces appearing in the 
block Numbers must appear in the order the 
surfaces are described on Card GQ, 

3CD(I) - flag to indicate end of input if blank, a comma 
means to continue. 

Continuation of Card GK in 12(15, Al) format is permissible. 

CARD GL (A4, A2, 1813) 

SI - hollerith SECT, 

DUM - dummy, 

IND(I) - the designation of each sector wliich describes 

the position of the sector relative to quadratic 
surfaces, 

+1: sector is on positive side of surface, 

-1: sector is on negative side of surface, 

0: surface is not needed to define sector. 

There must be a Card GL for each sector and references to 
quadratic surfaces must be in same order as they are listed 
on Card GQ. 

CARD GM (3A4, 10(15, Al)) (omit if NSTAT on Card GA = 2) 

NAM2 - hollerith REGI, 

DUM(2) - dummy, 

INPOO - a list of regions sector by sector in the block, 

BCD(I) - flag to indicate end of input if blank, a comma 
means to continue. 

Continuation with 12(15, Al) format is permissible. 





\ 


CARD GN (3A4, 

NAM2 

DUM(2) 

INP(I) 

BCD(I) 


10(15, Al)) (omit if NS TAT on Card GA = 2) , 

I 

- hollerith SURF, 

I \ 

dummy, 

same as for Card GK except for region input in¬ 
stead of material inpu^. 


CARD GO (A4, A2, 1813) (omit if NSTAT on Card GA = 2) 

SI - hollerith SECT, 

DUM -' dummy, 

IND(I) - same as for Card GL except for region input in' 

stead of material input. 

Repeats Cards GI to G0 for each block. 

-j-t- 


CARD GP (15, 16A4, A2) • 


N0BD 


total number of quadratic surfaces in ihe entire 
system, '■ 



hollerith characters ignored by the code. (Help¬ 
ful in identifying input at a later time.) 1 


CARD GQ (4 (DIO. 5> A4, IX, Al)) 

I 1 l 

C0F(J) - coefficient of the term, 

BCD1(J) - hollerith indicating which term of the equation. 

XJQ, YSQ, ZSQ, XZ, YX, YZ, XY, ZX, YZ, 

X, Y, Z, or blank are the possibilities,. 

BCD2(J) - a flag which indicates the quadratic equation 
continues. Any non-blank character ends ihe 
f|eld. The next function .must start on row card. 

Repeat Caras GQ until all surfaces have been described. 



of the input is shown in Table BUI. 


i 

i 


\ 
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TABLE BHI. (Cont’vJ.) 

FORTRAN 










B.3 Combinatorial Geometry Input Instructions 

The combinatorial geometry input data is read by the J0MIN 
subroutine, except for the region volumes VN0R(I), which are read by 
the GTVLIN subroutine whenever IV0PT = 3. For clarity of terminology, 
the terms ’’regions” and ’’media" have essentially the same meaning as 
in the 05R Geometry Package, but are constructed in a different manner. 
The term "zone” is the same as the "region” as defined in the original 
combinatorial geometry package. The term "body" has the same mean¬ 
ing as in the original combinatorial geometry package. 

CARD CGA (215, 10X, 10A6) 

IV0PT - option which defines the method by which region 
volumes are determined; if 

IV0PT = 0, volumes set equal to 1., 

IV0PT = 1, concentric sphere volumes are cal¬ 
culated 

IV0PT = 2, slab volumes (1-dim.) are calculated, 

IV0PT = 3, volumes are input by card, 

IDBG - if IDBG >0, subroutine PR is called to print re¬ 
sults of combinatorial geometry calculations 
during execution. Use only for debugging, 

JTY - alphanumeric title for geometry input (columns 

21-80). 


CARDS CGB (2X, A3, IX, 14, 6E10. 3) 

One set of CGB cards is required for each body and for the END 
card (see Table B IV). Leave columns 1-6 blank on all continu¬ 
ation cards. 


* 

Not operational. 
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ITYPE 

IALP 


FPD(I) 


specifies body type or END to terminate reading 
of body data (for example B0X, RPP, ARB, etc.). 
Leave blank for continuation cards, 

body number assigned by user (all input body 
numbers must form a sequence set beginning at 
1). If left blank, numbers are assigned sequen¬ 
tially. Either assign all or none of the numbers. 
Leave blank for continuation cards, 

real data required for the given body as shown 
in Table B IV. 


CARDS CGC (2X, A3, 15, 9(A2, 15)) 

Input zone specification cards. One set of cards required for 

each input zone, with input zone numbers being assigned sequen¬ 
tially. 

IALP - IALP must be a nonblank for the first card of 
each set of cards defining an input zone. If 
IALP is blank, this card is treated as a con¬ 
tinuation of the previous zone card, 

IALP = END denotes the end of zone description. 

NAZ - total number of zones that can be entered upon 

leaving any of the bodies defined for this input 
region (some zones may be counted more than 
once). Leave blank for continuation cards for a 
given zone. (If NAZ £ Oon the first card of the 
zone card set, then it is set to 5). This is used 
to allocate blank common, 

Alternate IIBIAS(I) and JTY(I) for all bodies de¬ 
fining this input zone. 

IIBIAS(I) - specify the "0R" operator if required for the 
JTY(I) body, 

JTY(I) - body number with the (+) or (-) sign as required 
for the zone description 
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CARDS CGD (1415) 


MRIZ(I) - MRIZ(I) is the region number in which the "I 

input zone is contained (I = 1, to the number of 
input zones). Region numbers must be sequen¬ 
tially defined from 1. 


CARDS CGE (1415) 

MMIZ(I) - MMIZ(I) is the medium number in which the 
,f lthtt i n p U t zone is contained (I = 1, to the 
number of input zones). Medium numbers must 
be sequentially defined from 1. 


CARDS CGF (7E10. 5) (omit if IV0PT 4 3) 

VN0R(I) - volume of the "I th '’ region (I = 1 to MXREG, the 
number of regions). 

Note: If ENDRUN is used to obtain collision density and track 
length per unit volume estimate of fluence, then a data 
statement in ENDRUN must give a relationship between 
region and media. In this case only one medium may be 
in a region. 
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TABLE BIV. (Cont'd.) 










B. 4 M0RSEC—Cross-Section Module Input Instructions 


CARD XA (20A4) 

Title card for cross sections. This title is also written on tape 
if a processed tape is written; therefore, it is suggested that the 
title be definitive. 


CARD XB (1315) 


Column 

Name 

Description 

5 

NGP 

the number of primary groups for which 
there are cross sections to be stored. 

Should be same as NMGP input in M0RSE, 

10 

NDS 

number of primary downscatters for NGP 
(usually NGP), 

15 

NGG 

number of secondary groups for which there 
are cross sections to be stored, 

20 

NDSG 

number of secondary downscatters for NGG 
(usually NGG), 

25 

INGP 

total number of groups for which cross 
sections are to be input. 

30 

ITBL 

table length, i. e., the number of cross 
sections for each group (usually equal to 
number of downscatters + number of up- 
scatters + 3), 

35 

ISGG 

location of within-group scattering cross 
sections (usually equal to number of up- 
scatters + 4), 

40 

NMED 

number of media for which cross sections 
are to be stored--should be same as MEDIA 
input in M0RSE, 

45 

NELEM 

number of elements for which cross sections 
are to be read, 
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CARD XB (Cont'd.) 
Column Name 


Description 


50 

NMDC 

number of mixing operations (elements 
times density operations) to be performed 
(must be ^ 1), 

55 

NC0EF 

number of coefficients for each element, 
including Pq, 

60 

NSCT 

number of discrete angles (usually 
NC0EF/2 Integral ), 

65 

ISTAT 

flag to store Legendre coefficients if 
greater than zero. 


CARD XC (1115) 


Column 

Name 

Description 

5 

IRDSG + 

switch to print the cross sections as they 
are read if > 0, if < 0 card sequence is not 
checked, 

10 

ISTR + 

switch to print cross sections as they are 
stored if > 0, 

15 

IFMU f 

switch to print intermediate results of p's 
calculation if > 0, 

20 

IM0M f 

switch to print moments of angular dis¬ 
tribution if > 0, 

25 

IPRIN f 

switch to print angles and probabilities if 
>0, 

30 

IPUN t 

switch to print results of bad Legendre co¬ 
efficients if > 0, 

35 

IDTF f 

switch to signal that input format is DTF-IV 
format if > 0; otherwise, ANISN format is 
assumed, 


Switches are ignored if IXTAPE < 0. 
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CARD XC (Cont’d.) 


Column 

Name 

Description 

40 

IXTAPE 

logical tape unit if binary cross-section 
tape, set equal to 0 if cross sections are 
from cards. If negative, then the pro¬ 
cessed cross sections and ether necessary 
data from a previous run will be read; in 
this case (IXTAPE < 0) no cross sections 
from cards and no mixing cards may be in¬ 
put. The absolute value of IXTAPE is the 
logical tape unit, 

45 

JXTAPE - 

logical tape unit of a processed cross- 
section tape to be written. This processed 
tape will contain the title card, the vari¬ 
ables from common L0CSIG and the per¬ 
tinent cross sections from blank common, 

50 

I06RT 

logical tape unit of a point cross-section 
tape in 06R format, 

55 

IGQPT 

last group (M0RSE multigroup structure) 
for which the 06R point cross sections are 
to be used NMGP). 


CARD XD (1415) (omit if IXTAPE £ 0) 

Element identifiers for cross-section-tape. If element identifiers 
are in same order as elements on tape, the efficiency of the code 
is increased due to fewer tape rewinds. 


CARD XE (omit if IXTAPE 4 0) 

Cross sections in ANISN format if IDTF * 0, otherwise, DTF-IV 
format. Cross sections for INGP groups with a table length ITBL 
for NELEM elements each with NC0EF coefficients. 
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CARDS XF (215, E10. 5) (NMIX cards are required, omit if 
' IXTAPE <0) 

KM - medium number (media numbers from 1 to 

MEDIA, see Card B, must appear on some XF 
card), 

KE - element number occurring in medium KM (neg¬ 

ative value indicates last mixing operation for 
that medium and at least one negative value is 
required for each medium), 

RH0 - density of element KE in medium KM in units of 

atoms/(barn cm). 


CARDS XG (15) (omit if I06RT * 0) 

NXPM - number of point cross-section sets per medium 
found on an 06R tape, 

= !, total cross section only, 

= 2, total and scattering cross section, 

= 3, total, scattering, and u *fission cross 
section. 

XCHEKR Card (415) Cross sections and cross-section input data may 
be checked independently of M0RSE utilizing XCHEKR. (5) The 
input to XCHEKR consists of the cross-section cards XA through 
XG preceded by a card as follows: 

IADJM - set greater than zero for an adjoint problem, 

MEDIA - number of cross-section media; should equal 

NMED on Card XB, 

NMGP - number of primary particle energy groups for 
which cross sections are to be stored; should 
equal NGP (or NGG if NGP = 0) on Card XB, 

NMTG - total number of energy groups for which cross 
section are to be stored. Should be equal to 
NGP + NGG on Card XB. 




B. 5 SAMB0 Analysis Input Instructions 


The following data are read from cards by SC0RIN: 


CARD AA (20A4) 

Title information—will be immediately output. 


CARD AB (815) 

ND 

NNE 

NE 

NT 


NA 

NRESP 

NEX 


NEXND 


number of detectors (set =1 if sO), 

number of primary particle energy bins to be 
used, (must be ^ NE), 

total number of energy bins (set =0 if £ 1), 

number of time bins for each detector (may be 
negative, in which case I NT j values are to be 
read and used for every detector) (set =0 if 
I NT I s i), 

number of angle bins (set =0 if £ 1), 

number of energy-dependent response functions 
to be used (set =1 if £ 0), 

number of extra arrays of size NMTG to be set 
aside (useful, for example, as a place to store 
an array of group-to-group transfer probabilities 
for estimator routines. If the subroutine 
ENDRUN which outputs fluence estimates from 
collision and track lengths is used, then this 
number must be at least MXREG + 2 (see Card I 
for MXREG), 

number of extra arrays of size ND to be set 
aside (useful, for example, as a place to store 
detector-dependent counters). 
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CARDS AC (3E10.4) (ND cards will be read) 

X, Y, Z - detector location. (If other than point detectors 
are desired, the point locations must still be in¬ 
put and can be combined with additional data 
built in to user routines to fully define each de¬ 
tector. ) 

Note that the distanc c between the above points and the XSTRT, 
YSTRT, ZSTRT values and the initial age, AGSTRT, will be 
used to define the lower limit of the first time bin. 


CARD AD (20A4) 

Title or units for total responses for all detectors. Will be used 
in columns 54 through 133 of the title for the print of these arrays. 

CARD AE (20A4) 

Title or units for each total response for all detectors. 


CARDS AF (7E10.4) 

Response function values. NMTG values will be read in each set 
of AF cards. Input order is from energy group 1 to NMTG (order 
of decreasing energy). 

Note: Cards AE and AF are read in the following order: AEj, 
AF},..., AF n , AE2 , AFj, ..., AF n , etc. NRESP sets of 
AE, AF cards will be read. 


CARD AG (20A4) (omit if NE s 1) 

Units for energy-dependent fluence for all detectors. 
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CARDS AH (1415) (omit if NE £ 1) \ 

t \ ! 

Energy group ^umbers defining lower limit of energy bins (in 
order of increasing group numbe?)~The NNE (if > 0) entry mi 
equal NGPQTN; the NE entry must be set to NMGP -i NGPQTG 
lor a combined problem,, or else NGPQTG or NGPQTN. 

\ 

CARD AI (20A4) (omit if j NT | ^ 1) ' 

I \ 

Units for time-dependent total responses for all detectors. ! 

CARD AJ (20A4) (omit if j NT j £ 1 or NE £ 1) ' 

Units for Line and energy-dependent fluence for all detectors. 

CARDS AK (7E10. 4) (omit if | NT | £ 1) 

NT values of upper limits of time bins for each detector (in ord 

aL 1 I!?f rt ' aS1,lg . tl ? e and detector number). • The values for each 
detector must start on a new card. | NT f values only are read 
if NT is negative. They are then used for every detector. 

CARD AL (20 A4) (omit if NA si) 

I i 

I i 

Units for angle- and energy-dependent fluence for all detectors. 

) 

CARD AM (7E10.4) (omit if NA £ 1) 

\ 

y* 1 ?.: 3 i u PP er Iim Ms of angle bins (actually cosine bins; 
the NA 11 * value must equal one). . 

Following the input for the SAMB0 analysis module, input cards 
for user-written routines INSC0R, S0URCE, and ENDKUN. 









\Input to PICTURE 


Card PA: Format (15) 


NUSE: 


The number of characters to reald on Card B to replace 
the standard values of ATABLE. Leave Card A blank 
and omit Card B if the standard ATABLE is desired. 
NUSE £ 50. 


Card PB: Format (50A1) (omit if NUSE = 0) 

ATABLE(I), 1=1, NUSE: The list of characters that are to be 
printed for each medium. For medium N, ATABLE 
(N+l) is printed. \ If N s 47, ATABLE (48) is printed. 
The standard values of ATABLE are: 


Medium Number 
0 (external void) 
1 through 9 
10 through 35 
36 through 46 

^ 47 (including 
internal voids) 


Character Printed 

1 through 9 
A through Z 

various special, 
characters 

(blank) 


GE0M input: Combinatorial geometry input. 
C£rd PC: Format (212, 18A4) 


= 0 After this picture, return to Card PC for another 
ICNT* picture with the same geometry. 

= 1 After this picture, return to Card PA to read in a 
ne^’ GE0M input. 

! = -1 Display the region geometry. 

= 0 Display the zone geometry. 

= 1 Display the material geometry. 

'TITLE(I), I = 1, 18: 72 characters to be printed as a title. 
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5.. 


Card PD: Format (6E10. 5) 



UL 


X, Y, and Z coordinates in the combinatorial geometry 
of the upper left corner of the picture. 


X 


LR 


LR 


LR 


X, Y, and Z coordinates in the combinatorial geometry 
of the lower right corner of the picture. 


NOTE: Card PD partially describes the plane of the slice by 

defining two points in the plane and designates the top, 
bottom, left and right sides of the picture. 


6. Card PE: Format (6E10. 5) 

u x ) 

I Direction numbers proportional to the direction cosines 
U Y / for the U axis of the picture. The U axis points down the 

TT 1 printed page in the direction the page moves through the 

U Z ) printer. 

v x } 

Vy \ Direction numbers for the V axis of the picture. The V 
( axis points to the right across the page. 

v z ) 


NOTE: Card PE completes the description of the plane of the 

slice by giving a line in the plane, also specifies the 
orientation of the picture on the output. 


7. Card PF: Format (215, 2E10. 5) 

NU: Number of intervals to print along the U axis (overrides 

DELU) . 

NV: Number of intervals to print clong the Y axis (overrides 

DELV) . 
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DELU: Spacing x ±n GE0M units) of intervals along the U axis. 
DELV: Spacing (in GE0M units) of intervals along the V axis. 

NOTE: All four entries are not required as input on Card PF, 

see below for explanation. 


Options: 

X LR= Y LR= Z LR=°- 

For this case NU and NV must be specified. In addition, either 
DELU or DELV must be specified. If the other is left blank, the 
code will produce an undistorted picture. If both DELU and DELV 
are specified the picture is likely to be distorted. standard 
printers give 10 characters to the inch across a line but only 6 
lines per inch down the page. Because of this DELV = . 6*DELU 
is necessary to produce an undistorted picture. 

X LR or Y LR or Z LR ^ °* 

If any one variable on Card PD is specified, the code will calculate 
the others to produce an undistorted picture. 

If both NU and DELU (or both NV and DELV) are specified, DELU 
(or DELV) will be ignored. 

The U and V axes may have arbitrary orientation. (If they are not 
orthogonal, the resulting picture will be distorted.) In Option 1, 
the first point will be at (X, Y, Z)ul> and the remaining points in 
the directions and at the distances specified. In Option 2, the 
range from Xul to X^r is ^-'vided into intervals and the calcu¬ 
lated points will be at the midpoints of the intervals. The first 
point will be 1/2 interval past (X, Y, Z)ul and the final point will 
be within 1/2 interval of (X, Y, Z)lr. If (X, Y, Z) L r does not 
lie on the U-V plane, or if the U and V axes are not orthogonal, 
the location of the final point is not readily predictable. 

The simplest method to obtain the correct results is to specify two 
diagonal corners of the plane of the slice on Card PD, with the top 
having the short dimension, if it is not square, and the side having 
the long dimension. Then, on Card PE, specify the U axis to be 
parallel to the edge of the slice with the large dimension (left or 
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right side), and the V axis to be parallel to the edge of the slice 
with the small dimension (top or bottom). Finally, let the only 
entry on Card PF be NV equal to maximum number of characters 
per line on your printer, this will provide the largest undistorted 
picture. 
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